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ABSTRACT 


The ability to make computer images more realistic is becoming more 
important as the hardware for producing such images is becoming less expensive 
and hence more available. The key to producing realistic images lies in the 
algorithms that can take full advantage of the hardware to produce them. In this 
study, we look at a prototype of a ray tracer, as presented in {Ref. 1]. Ray tracing, 
in combination with a global illumination model, currently provides the most 
realistic images that can be generated on general purpose computing hardware. 


The prototype was succesfully implemented on an IBM AT clone. 
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I. INTRODUCTION 


From the beginning of recorded history mankind has always had the need to 
create pictures. The reasons for creating these pictures ranges from the aesthetic, 
pretty pictures are nice to look at, to the functional, pictures can be an excellent 
way to communicate information. As mankind progressed, so did his ability to 
create pictures, although the techniques used to create pictures basically stayed 
the same. The advent of computers gave man yet another tool with which to 
create pictures. | 

The rapid increase in technology has made computer graphics a rapidly 
growing field. For the first time since man started drawing pictures, completely 
new techniques needed to be developed. In computer graphics the brush, paint 
and canvas are replaced by the mouse, algorithm and display. However even 
though the tools have changed, the same problems remain: how to make the 
picture look better, be it either more pleasing to the eye or to get the information 
across more clearly. 

Two of the most common and difficult problems in computer graphics are the 
hidden surface and lighting and shading problems. A large number of solutions 
exist to both of these problems. Very few solutions can be applied to both. One 


such solution is ray tracing. Ray tracing is the process of following an imaginary 


ray from a viewpoint through a pixel on a screen and into a scene to determine if 
it intersects any objects in the scene and then calculating the intensity of the pixel 
it went through based on the the final destination of the ray. As in most cases 
where one solution is found to several problems, that solution is seldom the best 
for all the problems it is applied to and so it becomes a matter of trade-offs. Such 
is the case with ray tracing. Among the hidden surface removal techniques, ray 
tracing is the least efficient being referred to as a brute force technique. In 
contrast, it has been labeled as one of the most elegant techniques in regards to 
lighting and shading (Ref. 2: p. 137]. Because of this latter fact, ray tracing has 
become an important technique in computer graphics. Ever since the idea behind 
ray tracing was suggested by Appel, numerous articles, studies, and 
implementations have been done on it. These in turn have spawned fruther 


extensions and modifications. [Ref. 1: p. 296] 


A. DEFINITION AND OVERVIEW 

The idea behind ray tracing lies in the theory that the light in our 
environment can be modeled as rays. After being emitted from a source, the rays 
are then reflected and refracted through a scene. Some of the rays eventually find 
their way to the eye where the scene is recreated (Figure l1.la). These light rays 
are emitted from light sources, such as the sun. An infinite number of light rays 
exist, but only a small percentage of them are received by us. To try and trace 


these rays from the source is computationally expensive. Appel suggested that 
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instead of tracing the rays from the source that they should be traced backwards 
from the viewer, thus dealing with only those rays that actually contribute to the 
scene (Figure 1.1b). [Ref. 1: p. 296] 

The basic ‘ai tracing algorithm is a very simple one and not difficult to 
implement. The basic algorithm is a hidden surface algorithm. All hidden surface . 
algorithms can be classified based on the coordinate system or space in which they 
are implemented. These are either in object space or in image space. The ray 
tracing algorithm falls under the category of image space. This category of 
algorithm is implemented in the screen coordinate system in which the objects are 
viewed. Unfortunately, the calculations are performed only to the precision of the 
scene representation, which generally provides poor resolution. The image space 
algorithms work by comparing every object in the scene with every pixel. Such an 
algorithm is computationally expensive. Ray tracing algorithms have three parts: 
a viewpoint, a raster screen, and a set of objects (Figure 1.2). In the algorithm, 
the viewpoint is along the positive z axis. From this point, a ray is shot into the 
scene through the center of every pixel on the raster. Each of these rays is then 
traced and compared against every object in the scene to determine if there is an 
intersection with any of them. It is in the determination of a possible intersection 
point that a ray tracer spends anywhere from 75 to 95 percent of its time. If 
there is an intersection, then the intensity at the pixel is determined using the 
intersected object’s attributes and an appropriate illumination model. If there is 


no intersection, then the pixel intensity is determined by the background 
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Figure 1.2 
Layout of Ray Tracing Scene [Ref. 1: p.297] 


13 


intensity. This procedure is then repeated for every pixel on the raster. When 
the ray tracer is used as a hidden surface algorithm, intersection testing stops 
after the first intersection. Extensions to the basic ray tracing algorithm, that 
showed its usefulness in implementing a global illumination model, were originally 
implemented by Whitted [Ref. 4] and Kay [Ref. 5 and 6]. In these extensions of 
ray tracing, additional rays are calculated, specifically the reflected and refracted 
rays, and then tested to see if they intersect with any objects in the scene. This 
process of generating new rays and tracing them to check for possible intersections 
is continued until the rays either leave the scene or stack space is exceeded. In 
such a case, the remaining rays are treated as if they had left the scene. This 
process, illustrated in Figure 1.3a, for a single ray with intersections is easily 
represented using the tree structure shown in Figure 1.3b. Here each node of the 
tree represents a ray surface intersection. At each node, at least one and 
sometimes two subbranches are generated. One branch of each of the reflected 


and refracted rays is generated from the point. [Ref. 1: pp. 190-296] 


B. ORGANIZATION 

This study is broken into three areas: data requirements, ray tracing 
methodology, and the intensity problem. The first section reviews the data needed 
for a lighting and shading modeler, hereafter referred to as a renderer, of which a 
ray tracer is an integral part. The second section reviews the actual process of 


tracing a ray through a scene to be rendered. The third section looks briefly at the 
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illumination problem and how it relates to the ray tracer. The concluding 
chapters present the implementation, and known limitations of the model along 


with areas of future research. 
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II. DATA REQUIREMENTS 


A. OVERVIEW OF THE DATA REQUIREMENTS 
The importance of the intersection routines in the ray tracer is apparent in 
the fact that a ray tracer spends 75 to 95 percent of its time determining 
intersections {Ref. 1: p. 297]. The key to determining intersections, however, lies 
in large part on the data used to describe the scene that is being rendered. 
Information is needed not only to describe the entire scene that is being rendered 
but more importantly to describe each object in the scene. Scene data is that 
information needed to completely describe a picture, i.e., the number, kind, shape, 
and color of any objects in the picture along with the background intensity and 
light source information. This information must be properly ordered and broken 
down. Falby [Ref. 3] suggested that a scene be broken into three categories: 
object, view, and light. Each of these areas is examined below in the context of a 
ray tracing algorithm. 
1. Object Data 
The data pertaining to each object in the scene can be grouped into two 
categories: polygon and bounding volume. The reason for this breakdown is 
twofold. First, each object in the scene is composed of polygons. They are the 


basic building blocks of the scene. Second, in order to reduce the number of 
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a 


intersection checks, it is necessary to set up a boundary around each object so 
that the ray tracer only performs the intersection checks in the regions that 
actually contain an object. Such a boundary is called a bounding volume. We 
examine the polygon data first. 
a. Polygon Data 

Since the main focus of the ray tracing algorithm lies in determining 
the intersections between the rays shot into a scene and the objects that make up 
a scene, and since each object is comprised of polygons, the problem is really one 
of determining the intersection points between the rays and the polygons. From 
the fundamentals of vector calculus, it is known that in order to determine the 
intersection between a ray and a polygon only the vertices of the polygon are 
needed as well as the direction of the ray and a point on the ray. Since the object 
is to be constructed of polygons, its vertices are known. Therefore, it is only 
necessary to ensure that these points are stored in some manner, such as a record, 
so as to be accessible to the ray tracer. In order to determine what the intensity 
of the pixel is through which the ray passes, it is essential that the characteristics 
of the object whose polygon was intersected be available. Since an object is made 
of polygons, they inherit the characteristics of the object. These characteristics 
also need to be readily accessible and, therefore, need to be stored in some 
manner. The following is a list of the basic object characteristics that need to be 
available: (1) the specular, diffuse, and transmission coefficients; (2) the Phong 


specular exponent; and (3) the index of refraction, see Table 2.1. [Ref. 3: p. 68] 
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TABLE 2.1: OBJECT DATA 


FIELD NAME VARIABLE NAME VALUE 


Polygon Vertices 

















Diffuse Coefficient (0-1) 
Specular Coefficient real (0-1) 
Transmission Coefficient real (0-1) 
Unit Surface Normal real (0-1) 
Phong Specular Exponent integer (0-200) 


Index of Refraction real 


b. Bounding Volume Data 

The major disadvantage of ray tracing is that it takes so much time. 
This is hard to avoid since it is so computationally expensive. It is essential, 
therefore, that more efficient techniques be developed to assist in reducing the 
number of calculations. Several techniques already exist with the bounding 
volume being the most effective [Ref. 1: p. 298]. In the description of ray tracing 
given so far, it has been stated that a ray is checked to see if it intersects with any 
object. Upon dissecting this statement further, a better understanding of the 
intersection problem can be realized. Unless some optimization is done, the ray 
tracing algorithm is forced to do the following. Each ray must be checked for a 
possible intersection with each object. Since each object is made up of polygons, 
then the ray must be checked for a possible intersection with each polygon. For a 
complicated object, such as a teapot. this requires a large number of checks and 
must be done for each object. The purpose for establishing the bounding volume 
lies in two facts. The first is that generally scenes are mostly background with 


just a few objects, hence very few of the rays actually hit anything. Therefore, 
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most of the intersection tests done are a waste of time. Second, a ray can only hit 
one object at a time. To have it process through the entire list of objects, when 
intersections with most of them can be eliminated, is needless. A bounding 
volume is, therefore, a method of enclosing each object in the scene in a simple 
containment vessel, which in effect creates a boundary around the object. Once 
this boundary is established, the number of overall intersection tests can be 
greatly reduced, as in the example of a teapot, which might easily have over a 
hundred polygons. If it is surrounded by a bounding box consisting of just six 
polygons, the number of intersection tests can be significantly reduced. In this 
situation, instead of having to test each ray against each polygon of the object, 
only those rays that penetrate the bounding volume need to be checked. Thus 
the bounding volume is a way to filter out unnecessary intersection tests by 
limiting the tests to those rays that are most likely to intersect an object. 

Just as the use of a bounding volume greatly increases the efficiency 
of the ray tracing algorithm, the use of the right kind of a bounding volume can 
improve upon that even more. In Rogers [Ref. 1], the bounding volumes 
suggested are a bounding box and a bounding sphere (Figure 2.1), each of which 
has advantages and disadvantages. The bounding sphere is much easier to 
implement although it is less efficient in reducing the target area than is the 
bounding box, see Figure 2.1. The bounding box, on the other hand, is 
computationally expensive to implement. The data needed to establish a 


bounding sphere is minimal. It only requires a center point for the object and a 
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Figure 2.la 
Obiect Surrounded by Bounding 


Figure 2.1b 
Object Surrounded by Bounding Box [Ref. 1: p. 
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radius that encompasses every point of the object. The bounding box, on the 
other hand, requires far more data in that the polygons that make up the box 
must be described. 
2. View Data 
In rendering any scene, estan information can be applied to the scene 
as a whole. This information is grouped together to form the view data. This 
data consists of the viewpoint position, a constant to prevent division by zero, a 
refraction index for the global medium, the ambient light intensity, the 
background color, and the scene dimensions, see Table 2.2. [Ref. 3: pp. 74-75] 
3. Light Data 
To support a lighting and shading model, it is necessary to include 
certain information on the light source for the scene. That information must 
include the position of the light source, its intensity, its type, geometry, and 


dimension, see Table 2.3. 


TABLE 2.2: VIEW DATA 


FIELD NAME VARIABLE NAME | VALUES 


| Viewpoint real 
| No Zero Constant real (<0) 
Global Refraction Index real 


Ambient Light real (0-1) 
Background Light Ib real (0-1) 
Scene Size integers 
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TABLE 2.3: LIGHT DATA 


F IELD NAME | VARIABLE NAME VALUE 


Light Position real 
Intensity real (0-1) 
Type enumerated (point distributed) 





Shape enumerated (circular, rectangular) 
Dimensions real 


B. DATA STRUCTURE FOR A RAY TRACER 
Falby [Ref. 3] suggested a data structure for a multi-illumination model 
renderer. That data structure, with minor variations, has been used in this study. 
In [Ref. 3] a complete derivation of the data is presented, so for the purposes of 
this study only a brief description is given here. Figure 2.2 illustrates the layout 
of the data structure as used in this study. This data structure essentially consists 
of arrays of records layed out in a hierarchical organization. Starting from the 
highest level it consists of the following: a picture record, an array of light records, 
an array of objects, an array of subobjects, an array of common part records, an 
array of polygons, and three arrays for the vertices. Each of these is now 
examined. 
1. Picture 
Picture is a single record which contains the view data mentioned earlier. 
The lights array is an array of records, with one record for each light 


source in the scene. Each record contains the light data mentioned above. 
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3. Objects 

The objects array is an array of records, with one record for each object 
in the scene. In this study, an object is the highest order item in a scene. Just as 
the scene is divided up into objects, each with its own bounding volume, so is 
each object broken down into subobjects, each with its own bounding volume. 

4. Subobjects 

The subobjects array is an array of records, with each array belonging to 
one object record. For example Figure 2.3a shows one object, a barbell, that is 
divided into three subobjects which are: the left weight, the right weight, and the 
bar. The record layout for this is as illustrated in Figure 2.3c. A subobject is the 
smallest item in the scene. Each object has at least one subobject. A subobject is 
composed of polygons or it is a sphere. Using Figure 2.3a as an example again, 
the left and right weights are spheres and the bar, instead of being a perfect 
cylinder, is composed of polygons and actually has an octagonal shape, Figure 
2.3b. Aside from containing a pointer to the common part record, examined next, 
and data for its bounding volume, it also contains information on the subobject 
type. This subobject type field indicates the geometry of the subobject, i1.e., it is 
either a sphere or a polygonal object, which is an object composed of polygons. 
This information is stored because different intersection routines are used for each 
object type. Currently a 1 indicates that planar intersection routines should be 


used and a O indicates spherical intersection routines should be used. 
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Figure 2.3a 
Subobjects [Ref. 7] 
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5. Common Parts 
A good illustration of common parts is found in an ordinary 
checkerboard, Figure 2.4a. In this figure one object exists--the checkerboard. This 
in turn has one subobject, itself. This subobject has two common parts: the white 
squares and the black squares. Table 2.1 listed the characteristics of an object 
and it is in the common parts record that these characteristics are stored. Each of 
these common parts records contains a pointer to an array of polygon records. It 
is through this arrangement that the polygons inherit the characteristics of the 
object. Therefore, the common parts array, also called the Cparts array, is an 
array of records, with one array per subobject, and each common parts record 
points to its own set of polygons, Figure 2.4b. 
6. Polygons 
The polygons array, too, is an array of records with one array for each 
subobject. This is the smallest physical item in the scene and the one against 
which the actual intersections are determined. 
7. Vertex Array 
The vertex array is an array of points that define the polygons that 


compose the subobject. 


This data structure as presented by Falby [Ref. 3] proved itself to be both 
flexible and easy to use. An example of a data base that used this structure and 


which was used in testing this ray tracer can be seen in Appendix B. 
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Figure 2.4a - Example of Object with Two Common Parts 
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Il. RAY TRACING INTERSECTION CONSIDERATIONS 


The methodology behind ray tracing is quite simple. However, it does require 
an understanding of the fundamentals of vector calculus and geometric optics. A 
discussion of the fundamentals is beyond the scope of this study. An overview of 


some of the fundamentals is in order. 


A. RAY TRACING MECHANICS 
By its very definition, ray tracing is simply the tracing, or following, of a ray 
from its source through space and determining any possible intersections that may 
occur between it and an object. The natural way to model a ray in order to do 
this tracing is by using vectors. A vector is not only a precise way to represent a 
ray but the basic operations on vectors in three space, addition, subtraction, dot 
product, and cross product provide the tools necessary to determine the 
intersections. These tools, along with other techniques found in vector calculus 
and geometric, optics provide the means to deal with the two problems 
encountered in ray tracing, i.e., the ray direction determination problem ans the 
intersection problem. 
1. The Ray Direction Problem 
Solving the ray direction problem is both the first and last step 
encountered in the ray tracing process. Determining the initial ray from the view 
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position, usually referred to as the view ray, is the simplest to solve. Every point 
in a coordinate system can be associated with a ray, and determining the direction 
of a ray between two points can be solved by using vector subtraction. The last 
step in the ray tracing process is determining what takes place when a ray 
intersects an object. This requires the application of the laws of geometric optics. 
Once a ray strikes an object, either one or two additional rays will be generated. 
These new rays are referred to as the reflected and refracted rays, Figure 3.1. The 
three basic laws of reflection and refraction are listed as [Ref. 8: pp: 32-33] 
1. The incident, reflected, and transmitted rays all reside in a plane, known 
as the plane of incidence, which is normal to the surface of the object. 
2. The angle of incidence is equal to the angle of reflection O, = ©. 


3. The incident and transmitted ray directions are related by Snells’ law: 


n.sin@. = n,sinQ.. 
t t t t 


An illustration of these laws is shown in Figure 3.1. Rogers |[Ref. 1: p. 367] 
provides a method for determining the direction of the reflected and refracted 


rays. The direction of r, the reflection ray, and p, the refraction ray are given as: 


where 
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Figure 3.1 
Creation of Reflected and Refracted Rays [Ref. 3: p.55] 
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Here k, is the ratio of refractive indices, ke is the Fresnel coefficient, v is the 
normal vector in the direction of the incoming ray, v’ is the unit normal vector in 
the direction of the incoming ray, n is the unit surface normal, and 7, and n, are 
the refraction coefficients for mediums the rays pass through. This is illustrated 
in Figure 3.2. 
2. The Intersection Problem 

As stated above the intersection computation is the most time 
consuming part of the ray tracing process. It is not that the process itself is so 
dificult but because several steps need to be done for each iteration. Two types 
of intersection computations are required to be performed: determining the 
intersection between a line and a sphere and determining the intersection between 
a line and a polygon. The first type is the simplest to solve and is why the sphere 
is generally used as the bounding volume. The calculation of the intersection 
point between a line and a sphere involves solving the equation for the line and 


the sphere simultaneously. The sphere is defined by the equation 
(rz - a) + (y- f) + (2-8 =r (3.1) 


where (a, @, 6) is the center point, r is the radius, and (z, y, z) is a point on the 


sphere. The line is defined by the parametric equations 
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Figure 3.2 
Direction of Reflected and Refracted Rays [Ref. 3: p.57] 
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r= of bo 1) = bt UC eeer (3.2) 


where (Z,, Yg; 29) is a known point on the line, and a, b, and ¢ are coefficients 
from the vector at + b7 + ck which is parallel to the line. These equations must 
first be solved for t. The solution to this provides two answers. First, it indicates 
whether or not an intersection actually takes place. Second, if there is one, it 
indicates how many intersections, either one, in the case where the line is tangent 
to the sphere, or two, in the case where it actually enters the sphere. In the case 
of two intersection points, a check must be done to determine which is closer to 
the origin of the ray. 

The intersection between a line and a polygon is more involved. This 
problem is comprised of two parts: determining the intersection point between a 
line and a plane, and determining whether or not the intersection point lies within 
the polygon. Like the line-sphere intersection problem, this one also involves the 
solving of two equations simultaneously. The first of these is the equation for a 


plane which is defined as 
Az + By + Cz=D (3.3) 


where A, B, C, and D are constants and (z, y, z) is a point on the plane. The 
second equation used is the parametric equation that defines a line, Eq. 3.2. The 
solution of this requires first substituting the equations for the line into the 
equation of the plane. The result is an equation in t, which when solved and 


substituted back into the equations for the line, provides the intersection point 
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between the line and the plane. Once this point is calculated, it is then necessary 
to determine whether or not it lies within the polygon. A general way of solving 
this is to simply determine the relationship between the intersection point and 
each edge of the polygon. The point that lies on the inside of each edge also lies 
within the polygon. If it fails the test for any edge, then it lies outside the 
polygon. The drawback to this approach, is that it only works for convex 


polygons. In this study we assume all polygons are convex. [Ref. 9]. 


B. THE RAY DATA STRUCTURE 

Rogers |Ref. 1] suggests a data structure for a ray in a ray tracer. It is that 
data structure which is used in this study. Table 3.1 lists the data used to model 
each ray. We examine each item of this structure as adapted from Rogers [Ref. 1: 


p. 373]. 


TABLE 3.1- RAY DATA 


FIELD NAME VARIABLE NAME VALUE 






































Ray Type type enumerated or coded 
Ray Origin Te Yee real 

Ray Vector eee real 

Source Ray Type Stype enumerated or coded 
Intersection Flag flag boolean or coded 
Object Index obj idx integer 
Subobject Index | subobj idx integer 


Common Part Index | Cpart idx | integer | 



























Polygon Index polygon idx integer 
Intersection Point Ses real 
Distance real 
Transmitted Intensity real 





Specular Intensity real 


395 


1. Ray Type 
The ray type field identifies a ray as either a view ray, a reflected ray, or 
a refracted ray. The values put in this field are generally of an enumerated type 
and consist of reflected, refracted, view, and none. 
2. Ray Origin 
The ray origin field contains the point that identifies the position from 
which the ray originated. For instance, if it is the view ray, its point of origin is 
the view position. If it is a reflected or refracted ray, its origin is the intersection 
point that it originated from. 
3. Ray Vector 
The ray vector field contains the vector heading of the ray. 
4. Source Ray Type 
The source ray type field contains the ray type of the source ray for this 
particular ray. For instance, a view ray does not have any source ray as it is the 
starting ray for the process. Hence, none is in the type field. If the view ray 
intersects an object and both a reflected and refracted ray are generated, then the 
source ray for both of them is the view ray. Likewise if the reflected or refracted 
ray hits an object and creates further rays then it becomes the source ray for 
those rays it creates. 
5. Intersection Flag 
Originally the intersection flag is set to false and it is only set to true 


when there is an intersection between this ray and an object. 
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6. Object Index 
The object index provides an index into the array of object records 
making it possible to select any object easily. 
7.  Subobject Index 
The subobject index provides an index into the array of subobject 
records and helps uniquely identify each subobject. 
8. Common Part Index 
The common part index provides the index into the array of common 
part records uniquely identifying each common part record. 
9. Polygon Index 
The polygon index provides the index into the array of polygon records 
and uniquely identifies each polygon. 
10. Intersection Point 
The intersection point field holds the position of the intersection point 
between the current ray and an object. 
11. Distance 
The distance field contains the distance between the current ray’s point 
of origin and its point of intersection. 
12. Transmitted Intensity 
The transmitted intensity field contains the red, green, and blue 
intensity values, in a range between 0 and 1, of the light that is incoming along 


any refracted ray that this ray produces. 
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13. Specular Intensity 
The specular intensity field contains the red, green, and blue intensity 
values, in a range between O and 1, of the light that is incoming along any 


specular ray that this ray produces. 


C. INTERSECTION METHODOLOGY 
1. Intersecting a Planar Polygon 
a. Generating the Initial Ray 
The generation of the initial or view ray is shown in Figure 3.3. 
This ray, p — v, 1s calculated using vector subtraction. The two vectors used are 
the ones associated with the points for the view position, v, and the pixel through 
which the ray passes, p. The ray associated with the view position is to be 
subtracted from the ray associated with the pixel. 
b. Intersecting the Bounding Volumes 
After each ray is generated, each object in the object list is checked, 
one at a time, to determine whether or not the ray strikes any of the bounding 
containers of the objects in the scene. In our implementation, the bounding 
container is the sphere, which is reduced to a bounding circle, C, see Figure 3.4. 
his circle’s radius is the same as the radius of the bounding sphere, S, and it lies 
on the plane, P, which contains the center point, g, of the bounding sphere. The 
inverse of the incoming ray, ¥, is the surface normal of this plane. It is not 


necessary to determine where on the bounding sphere a ray hits since at this stage 
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we are only interested in just finding out whether or not it hits it. Because we are 
not concerned with where the ray strikes the bounding sphere, a bounding circle is 
used. Determining the intersection with a circle requires less work than 
determining the intersection with a sphere. 

The first step in tracing the ray begins by taking each object and 
constructing a bounding circle. Once this circle is constructed, the intersection 
point of the view ray with it is calculated. Then the distance between the 
intersection point and the center point is determined and compared to the radius 
of the bounding circle to see if it falls within the circle. If it does fall within the 
circle, this indicates that the ray intersected the bounding volume. This process 
Hien needs to be repeated for each object. If an object’s bounding circle is hit, 
this process must then be repeated for each subobject of that object. 

c. Intersecting the Polygon 

Once a particular subobject is identified by the bounding volume 
tests, the common parts list is processed, one record at a time. Each one of these 
common parts records contains a pointer to those polygons that make up the 
object being rendered. This list of polygon records is then processed after the 
bounding volume processing is completed. The entire list needs only to be 
processed until an intersection is found. The processing involved in this is the 
most computationally expensive part of the entire ray tracing algorithm. This 
computation consists of three steps: determining the orientation of the polygon, 


calculating the intersection point between the ray and the plane that contains the 
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polygon, and determining whether or not the intersection point lies on the 
polygon. 

(1) Establishing the Orientation. First, each polygon needs to be 
checked for a correct orientation. This is a straightforward step carried out by . 
calculating the angle between the surface normal of that polygon and the inverse 
of the view ray. If it is 90 degrees or greater, it is facing in the wrong direction to 
be intersected. If it is less than 90 degrees, the next step is to determine whether 
the ray intersects the particular plane that that particular polygon lies on. The 
first step in doing this is to determine the equation for the plane in which that 
polygon lies. 

(2) Intersecting a Plane. If the correct orientation exists for a 
polygon to be intersected, the next step is to define the plane containing the 


polygon of interest. The equation for a plane was given earlier as 
Ar + By + Cz=D 


where A, B, C, and D are constants and can be calculated by the following 


equation. 


D = —2,(Y_23 — Y32_) — 2q (¥32, — ¥y23) — 23 (Yy2_ — Y2?4) 
where (z,, ¥,, 2), (Zo: Yq: 22), and (23, y3, Z,) are points of the vertices of the 
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polygon. Once the equation for the plane is known, it must then be solved 
simultaneously with the equation of the line representing the ray. 

(3) Location of a point with respect to a polygon. Once the 
intersection point is determined, it needs to be checked to see whether or not it 
lies within the polygon. This is the most computationally expensive part of the 
process. This process, see Figure 3.5, requires that a plane, P, called the bounding 
plane, containing two vertices, for example A and B, of an edge and another 
arbitrary point, A, not on the polygon or its plane, be constructed. This must be 
done for each edge. Once the bounding plane is constructed, the point in question 
must then be checked to see whether or not it lies to the polygon’s side of the 
plane that now contains the edge of the polygon. This is done by plugging in 
another vertex of the polygon into the equation for the plane that was just 
constructed, and then plugging in the intersection point. If the results from these 
two equations have the same relationship, i.e. if the sign is the same, then the 
intersection point lies on the polygon side of the bounding plane. This check 
must then be repeated for each edge of the polygon. In order for the intersection 
point to lie on the face itself, it must be found to lie on the polygon’s side of each 
edge. At this point, it is determined whether or not the ray strikes the object. If 
the ray does strike the object, the pointers to this polygon must be stored in the 
ray structure along with the intersection point. This information is needed later 
on when determining the direction of reflected and refracted rays as well as the 


intensity of the light that is reflected from this position. 
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Location of a point with Respect to a Polygon [Ref. 9] 
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2. Intersection of a Sphere 

The sphere is the easiest object to work with in a ray tracer. Since the 
sphere can act as its own bounding volume, the center point and radius are 
already available in the subobject record. This eliminates the need for the polygon 
array. The center point and radius are the only information necessary to model a 
sphere for a ray tracer. Determining the intersection of a line and a sphere is 
nothing more than the simultaneous solving of their equations for the variable t . 
The solution to this gives a quadratic equation in t which can then easily be 


solved. 
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IV. THE INTENSITY PROBLEM 


One of the strong points of the ray tracing isles is that a global 
illumination model can easily be integrated into it. In fact, ray tracing and global 
illumintion models seem to naturally complement each other. A_ global 
illumination model takes into account all of the light sources in the scene in 
calculating the intensity at each point. This means taking into account the 
ambient light that exists in a scene, light that comes directly from a particular 
light source(s), and light that is reflected off an object(s). It also includes the 
coefficients necessary to model the way an object reacts with light. A great deal 
of work has been done in this area. The most notable model is the Whitted 
illumination model, and it is the one that has been implemented here [Ref. 1: pp. 
365-366]. The Whitted algorithm is based on the three models shown in Figures 


4.1, 4.2, and 4.3. These models will now be examined more closely. * 


A. LOCAL ILLUMINATION MODELS 
1. Diffuse Reflection Model 
The first of these models is a perfect diffuser. Such a model is governed 
by Lambert’s cosine law. This law states that the intensity of light reflected from 
a perfect diffuser is proportional to the cosine of the angle between the light 


*The contents of this chapter are close adaptations from Rogers (Ref. 1] and Falby [Ref. 3]. 
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direction and the normal to the surface. This can be expressed mathematically as 


I = Lk ,cosO 0<O<— 
where J is the reflected intensity, J, is the incident intensity trom a point light 
source, k, is the diffuse reflection constant, unique to each object, and © is the 
angle between the light direction and the surface normal see Figure 4.1. Since the 
diffuse reflection coefficient &, varies from material to material and is also a 
function of the wavelength of the light, it is often easier to just assume it a 
constant for simple illumination models. [Ref. 2: p. 312] 
2. Specular Reflection Model 

The second model illustrates the characteristics of specular reflection 
which is directional, unlike diffuse reflection. This means that the greatest 
intensity of the specularly reflected light can only be seen if the view angle 
coincides with the reflection angle, Figure 4.2. The further off the viewing angle is 
from the reflection angle, the dimmer the intensity becomes. Because of the 
complex physical characteristics of specularly reflected light, an empirical model 


due to Bui-Tuong Phong is usually used for simple illumination models |Ref. 10]. 


This is expressed mathematically as 


where w(O., A), the reflection curve, gives the ratio of the specularly reflected 
hght to the incident light as a function of the incidence angle ©. and the 


wavelength A. Because w(O,, A) is such a complex function, it is frequently 
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Figure 4.3 
Global Illumination Model 
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replaced by an aesthetically or experimentally determined constant k, which then 


yields 
i = [k.cos" a 


Also, n is a power that approximates the spatial distribution of the specularly 
reflected light. Typically a value of 200 for n is used to model a very shiny 
surface and a value of 10 is used for a dull surface [Ref. 3: p. 72]. [Ref. 1: pp. 313- 
315] 
3. Combined Model 

If just point sources are assumed, as in the two models just discussed, 
any object not receiving light directly from the source appears black. In order to 
properly render a scene, it is also necessary to take into account ambient light, the 
light that is reflected off other surfaces. Including a model for ambient light into 
the intensity calculations is not feasible. Ambient light represents a distributed 
light source and as such is a very complex function. Therefore it is treated as a 
constant diffuse term and linearly combined with the other terms. Also not 
included in the above model is the effect that distance has on light. It is well 


known that the farther away an object or light source is, the dimmer it gets. The 


actual formula to produce that proper attenuation is ——, where the intensity of 
aa 


light decreases as the square of the distance from the source increases. However it 


has been shown that linear attenuation can actually produce more realistic results. 


With these two additions the complete model now looks like: 
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if 





(k,cosO + k.cos a) 
ductile 


where J, is the incident ambient light intensity, &, is the ambient diffuse 
reflection constant and A is an arbitrary constant that can be used to adjust the 
results. [Ref. 1: p. 313] 

The complete model just presented can now be modified to better fit in 
with a ray tracer. Recalling the formula for the dot product of two vectors allows 


writing the cos© as ri-Z and writing cosa as R-S which gives us: 


I 


d+kK 





[kK (aL) + k,(R-S)"] 


So far we have only examined the case where just one light source is present. If 
there are several light sources, the effects are added linearly, and the equation now 


becomes: 


{. 


d 





I=1k +3 


m OL 
aI 


a 7 eal 2) + k,(R,S)"] 


This then is the complete local illumination model. [Ref. 1: pp. 312-316] 
B. GLOBAL ILLUMINATION MODEL 
The complete local illumination model just presented forms the basis tor the 


algorithm that was implemented for this study [Ref. 1: pp. 363-378] see Figure 


4.3. 
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T=k1I,+, EI, ne i + kD, 5 R | =p Mee ted 
j i) 
In the above equation k,,k,,k,, and k, are the ambient, diffuse, specular 


reflection, and transmission coefficients, all of which have values between 0 and 1. 


I, I,, I,, and I, are the intensities of the ambient light, the specularly reflected 
3 


light, the transmitted light, and the light directly from a light source. These also 
hold values between O and 1. The remaining variables ni, L FP S, and R, are the 
surface normal at the intersection point, the direction of the jth light source, the 
local sight vector, and the local reflection vector from the jth light source. A 
careful comparison between this model and the complete local illumination model 
reveals that the only new terms are the J,and J, terms. These are the terms used 
to account for the light that comes in along the reflected and refracted rays that 
originated at this point. J, holds the intensity for the reflected ray and J, holds 
the intensity for the refracted ray. These two values in turn are calculated using 
the exact same model. For the last intersection point in the scene, the one whose 
reflected and refracted rays do not intersect anything, J, and J, are set to 0. The 
k, and k, terms are coefficients included to better model the way this object reacts 
with the light incoming along the reflected and refracted rays. 

This then is the complete global illumination model used in this study. It is 
the simplicity of this algorithm that makes it so easy to understand and 


implement. In essence, it is saying that the output intensity is nothing more than 


a0 


a sum of all possible light sources with the coefficients determining the intensity of 


light that comes from a particular object. 
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V. RAY TRACING ALGORITHM 


Rogers suggested an algorithm for ray tracing [Ref. 1: pp. 374-377]. It is that 
algorithm that has been the basis for this study. The following is a description of 


that algorithm as it has been implemented here. 


A. TRACING THE RAYS 

To begin the ray tracing process, the first thing done is the determination of 
the direction of the view ray.* The ray data, mentioned in chapter III, is then 
initialized. After the view ray is generated and the ray data initialized, the ray, 
which is represented by this ray data, is pushed onto the stack, which is used to 
implement the ray tracing tree. The process then moves into the actual ray 
tracing loop. 

Once in the loop, the stack is checked to see if it is empty. If it is not empty, 
the stack is then popped to access the ray information. The first thing checked is 
whether or not the intersection flag is set. If it is, that indicates that that 
particular ray has already been terminated (by hitting an object), and the process 


of determining the intensities begins. 


*Each ray, is modeled as a vector, and is converted to a unit vector immediately after its 
determination. 
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If the intersection flag is not set, a new ray has been generated, either a new 
view ray or a new reflected or refracted ray. All of these rays are grouped under 
the more general title of a shooting ray. At this stage, the ray must be sent 
through the intersection procedures to determine whether or not it intersects any 
object. The intersection routines start at the highest level of the picture and step 
through the linked list of objects, subobjects, and common parts to the actual 
polygons. 

At the top level of the intersection routines, each one of the objects is 
checked. First, they are checked to see if there is an intersection with the object’s 
bounding circle. If there is an intersection, the distance between the intersection 
point and the origin of the ray is calculated and sllacad in the ray data. Second, a 
check is done to insure that the objects lie in front of the ray’s origin. In front 
refers to those objects that lie in the direction of the shooting ray from its point of 
origin. Since the ray is being modeled by a vector, which only indicates direction, 
every object along the line described by the vector and the origin of the ray is 
considered, see Figure 5.1. The way to test for this is to generate a test vector 
between the origin of the ray and any intersection points of the shooting ray, 
eliminating the intersection points that lie in the opposite direction from the 
shooting ray. 

As each object is processed, the objects whose bounding volumes are 
intersected by the shooting ray are processed further to see if any of their 


subobjects are also intersected. The same basic procedure for finding an 
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intersection is used here as is used for finding an intersection for the objects. If the 
object and the subobject are one and the same, this check need not be done. 

For each subobject whose bounding volume is intersected, the next step is to 
check the type of geometric object, i.e., a sphere or a polygonal object. This needs 
to be done in order to determine which intersection routines are needed. In the 
case of a sphere, it becomes just a matter of solving two simultaneous equations, 
discussed earlier--one for a line, the other for the sphere. When dealing with a 
polygonal object, it becomes more complicated because each face of the polygon 
must be checked. The first thing that needs to be done is to check the orientation 
of the face in question. When the polygon does have the right orientation, the 
intersection process continues. The first step is to determine the equation of the 
plane that contains the polygon. This plane is calculated from any three vertices 
of the face in question. Once the plane has been established, the intersection point 
between it and the shooting ray is calculated. With the intersection point thus 
established, the next step is to determine whether or not the point lies on the face 
of the polygon. If the intersection point is found to lie on the face of the polygon, 
then this point is saved and placed in the ray data. 

If no intersection is found and the shooting ray is either a reflecting or 
retracting ray, then it is discarded and the intensity calculation process begins. if 
the shooting ray is a view ray, then the intensity is set to the background 


intensity which is then displayed and the next ray is generated. 
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When there is an intersection, the first thing to be checked is whether or not 
there is enough room on the stack. Since the stack holds only part of the ray tree 
at any one time, it need only be long enough to contain the longest anticipated 
branch. A particular branch of the ray tree is terminated when both the reflected 
and Peracted rays at an object intersection leave the scene or when the available 
stack length is exceeded. When both rays leave the scene, their contribution to 
the illumination at the source ray is zero. When the available stack length is 
exceeded, the algorithm calculates the illumination at the source ray using only 
the ambient, diffuse, and specular reflection components at the source ray 
intersection. An extension algorithm is given in Rogers whereby the algorithm can 
be extended one additional depth in the tree without exceeding the maximum 
stack depth. However, the implementation of that was not necessary here. When 
the stack does get full, it becomes a matter of calculating the intensity at that 
intersection point and setting the amperes value J, or J, in the source ray. 
(Ref. 1: p. 372| 

When the stack is not full, then the distance between the source point of the 
shooting ray and the intersection point is determined and placed in the ray data. 
The ray is then placed back on the stack. Once that is done, any reflecting and/or 
refracting rays emanating from the intersection point are determined, their ray 
data initialized, and then placed on the stack as new rays with the reflecting ray 
being placed on first. It is important to keep this order of rays in mind because it 


is necessary to know the number of rays to pop when setting the intensities of the 
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source rays. With these new rays in place on the stack, the program returns to 
the beginning of the ray tracing cycle. In the absence of reflecting or refracting 
rays, then the first ray popped is the view ray. Now since this ray already has its 
intersection flag set, the intensity at the point of its intersection is calculated and 
displayed. If the ray popped is a reflective or refractive ray, then it becomes the 
new shooting ray. This cycle continues until either no more reflecting or 
refracting rays are produced or until the stack becomes full. This process is 


summarized by the pseudocode description in Figure 5.2. 


B. DETERMINING THE INTENSITY 

If a ray’s intersection flag is found set at the beginning of the ray tracing 
process, that ray is sent into the intensity processing part of the ray tracer. The 
first step in this process is to take the ray data and to determine the intensity at 
that ray’s intersection point. This intensity is then divided by the distance 
between the ray’s point of origin and its intersection point in order to properly 
attenuate it. This process is done for each of the primary colors -- red, green, and 
blue. If a view ray is being considered, that means that it was the last ray on the 
stack. Therefore, it is the final intensity to be calculated and hence is the actual 
one displaved. If it is a refracted ray. then the intensity just caiculated becomes 
the J, intensity in the source ray for the ray currently being examined. If it is a 
reflected ray, then the intensity just calculated becomes the J, intensity in the 


source ray for the ray currently being examined. The stack is set up so that the 
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READ IN DATA FILE 


for Y := 1 to MAX ROWS do 
for X := 1 to MAX COLUMNS do 
INITIALIZE VIEW RAY 
PUSH RAY ONTO STACK 
repeat 
POP RAY FROM STACK 
if INTERSECTION FLAG SET then 
CALCULATE INTENSITY 
else 
CHECK FOR INTERSECTION 
if INTERSECTION FLAG SET then 
if STACK EXCEEDED then 
CALCULATE INTENSITY 
else 
PUSH RAY BACK ON STACK 
CALCULATE REFLECTED RAY 
CALCULATE REFRACTED RAY 
if REFLECTED RAY EXISTS then 
INITIALIZE REFLECTED RAY 
PUSH REFLECTED RAY ON STACK 
end if 
if REFRACTED RAY EXISTS then 
INITIALIZE REFRACTED RAY 
PUSH REFRACTED RAY ON STACK 
end if 
end else 
else 
if (CURRENT RAY TYPE = VIEW RAY) then 
SET INTENSITY TO BACKGROUND COLOR 
end else 
until STACK EMPTY 
DISPLAY Bie 
end FOR-X 
end FOR-Y 


Figure 5.2 - Pseudocode Description of the Ray Tracing Process 
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source ray is on the bottom of the ray tuple with the reflected ray above it and 
the refracted ray above that. In order to set these values, it is necessary to pop 
the stack to gain access to the source ray. Once these values are set in the source 
ray, the program returns to the beginning of the loop and pops the next ray off 
the stack discarding the last ray as it is no longer needed. (Ref. 1: p. 377] 

The intensity algorithm [Ref. 1: p. 377] although very simple in design, 
depending on the number of light sources in the picture, can also become a time 
intensive part of the ray tracing program. The entire ray data set, listed in Table 
3.1, is sent into this routine as well as the pointers to the object and light source 
list. As in the intersection routine, rays are generated here. In Rogers, they are 
referred to as shadow feelers and the same term is used in this study. These 
shadow feelers are those rays represented as the vectors from the point of 
intersection to the light source, see Figure 5.3. They are used to determine the 
intensity contributed to that point from that light source. Once these rays are 
generated, they also pass through the intersection routines in order to determine 
which objects, if any, the light rays pass through en route to the intersection 
point. The first test that must be done is to determine if any of the objects passed 
through are opaque. If anv are opaque. then no light reaches the intersection 
point irom chat light source. That point is then considered to be lying in deep 
shadow with respect to that light source. If none of the objects intersected by the 
light ray is opaque, then the light intensity needs to be attenuated accordingly for 


each of them. This attenuation entails multiplying the intensity at each point by 
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Shadow Feelers 





Figure 5.3 
Shadow Feelers [Ref. 3: p. 64] 
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the transmission coefficient of the object. This process then needs to be repeated 
for each light source. A running total of the intensities is maintained to be 
included in the final calculation. The calculation takes into account the ambient 
light and the light that comes in along the reflection and refraction rays. This 
process then produces the final intensity. If the input ray type is a view ray then 
it is displayed. If it is a reflected ray, it becomes the J, value in the source ray. If 


it is a refracted ray, it becomes the J, value in the source ray. [Ref. 1: pp. 376-377] 
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VI. IMPLEMENTATION 


‘The prototype was written in BORLAND’s Turbo Pascal and implemented 
on an IBM AT clone. The program is 2500 lines long and takes two hours to 
generate a scene of 200x200 pixels. The scenes generated on the AT were then 
displayed using the RGB monitor of a Silicon Graphics IRIS 2400 graphics 
workstation. 

The main focus of this study was to develop a prototype ray tracer, which by 
itself is just a hidden surface removal technique. The secondary consideration was 
to integrate an illumination model into the ray tracer. Because of this focus, 
more time was spent examining the ray tracing algorithm than any of the global 
illumination models that could have been integrated along with it. 

The top three scenes in Figure 6.1 tested the ability of the prototype to 
perform as a hidden surface remover. The program proved successful in this area. 
For these scenes, a stub was used in place of the illumination model and each 
scene was lit using only ambient light. From left to right the scenes show: An 
unobstructed view of the three objects, described in a later section: the cube and 
sphere suspended above the floor but with the cube partially blocking the sphere; 


the cube and the sphere sunk part way into the floor with the cube still in front. * 


*These were the only scenes generated to test for hidden surface removal. The remainder of 
the tests were done trying to integrate a global illumination model. 
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The testing of the global illumination model integrated into the ray tracer 
unveiled some problems. The two bottom scenes in Figure 6.1 are representative 
of the successful results. The first problem discovered was that attenuating the 
intensity by the full distance between the origin of the view ray and its 
intersection point produced totally blackened objects. The results shown in 
Figure 6.1, the bottom row, were obtained by either dividing the intensity by two, 
see the scene on the left, or by not dividing it at all, see the scene on the bottom 
right of the object. The second problem can be clearly seen by the black line that 
runs up the center of the floor on the bottom right scene. This result along with 
thosé test that generated shadows, not shown, indicated that the intensities for 
the floor were reversed. The black line, clearly seen in color Figure 5 and just 
vaguely visible in color Figure 4, is actually specular reflection and should be 
much brighter than the rest of the floor. In those scenes where shadows were 
generated the shadows were also brighter than the rest of the floor--just the 


reverse of what it should have been. 


A. INEVe 

The test data used in this study produced the scenes shown in Figure 6.1. 
This test data was in the form of a sequential dle with the data structure outlined 
in Figure 2.2. The data had in it one picture record, one light in the lights array, 
and three objects in the objects array. The first object is a cube which contains 


one subobject. This subobject has one common part. The common parts record 
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then points to the six polygon records that where used to construct the cube. 
Each one of these records has its own vertex arrays. The second object is a 
sphere. This object has one subobject, and the subobject has one common part. 
Since a sphere can be its own bounding volume and since it is not constructed of 
- polygons, then the object chain for the sphere need go no further than the 
common parts array. The third object is the floor.* This consists of one 
subobject, one common part and one polygon. All objects in the scene are opaque 


and have a highly reflective surface. 


fm OUTPUT 

The output generated by the ray tracer is in the form of a bitmap, with 
values for each of the red, green, and blue components. These values range from 
Oto 1. To display these on the RGB monitor of the IRIS, each red, green, and 
blue component is then multiplied by 255 and assigned an index in the color 


table. 


*The floor is at a 10 degree angle to the screen to provide a better perspective. 
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VII. CONCLUSIONS 


A. AREAS OF FUTURE RESEARCH 

The ray tracer is a powerful tool in computer graphics. In its original design, 
it produced the finest rendered pictures at that time. Since then there have been 
numerous extensions, among the most widely known are the ones by Phong, Blinn 
and Newell, Kay, and Whitted each of which have further enhanced the 
performance of the ray tracer [Ref. 2: pp. 343-344]. There are two main areas for 
future research: global illumination models, and intersection algorithms. Both of 
the areas are of great interest in the graphics world. Since ray tracers and global 
illumination models can be integrated so easily, working on either problem would 
undoubtedly lead to insight into the other. This would also be very easy to do 
because ray tracing naturally lends itself to a modular design making it easy to 
establish hooks for the testing of a large number of algorithms, both illumination 


and intersection. 


B. CONCLUSIONS 

We have examined the three major areas of ray tracing: the scene data 
needed, the intersection problem, and the intensity problem. The data structure 
used was adapted from [Ref. 3] and proved to be useful. The intersection 
problem, although involved, is not complex. The algorithms used in this 
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implementation are simple. The inclusion of a simple global illumination model, 
was easy to integrate and provided fair results. The ray tracer provides an 
excellent test bed program and implemented can provide a useful tool to study 
numerous problems, not only in lighting and shading but also in intersection 


determination. 
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APPENDIX A — SOURCE LISTINGS 


DECLARATIONS 


const 

{ THIS IS THE CONSTANTS DECLARATION SECTION } 
maximum size of stack = 100; 
initial pixel _z = 0; 


type 


{ TYPES DECLARATION AREA } 
vertices array = array [1..4] of real; 


polygon rec = record 
num vertices : integer; 


vertice xX, vertice y, vertice 2 : vertices array; 
surface normal x, surface normal y, surface normal 2: real 
end; 


polygon ptr = “polygon array; 
polygon array = array [1..6] of polygon rec; 


common part rec = record 
K ar, K ag, K ab: real; { AMBIENT DIFFUSE COEFFICIENT } 
K dr, K dg, K db:real; { DIRECT DIFFUSE COEFFICIENT } 
K sr, K_sg, K sb: real; { SPECULAR COEFFICIENT } 
K tr, K tg, K tb: real; { TRANSMISSION COEFFICIENT } 
obj refraction index : real; { OBJECTS REFRACTION COEFFICIENT } 
obj phong exp’ : integer, { PHONG’S SPECULAR EXPONENT } 
num polygons : integer; 
polygons : polygon ptr 
end; 


common part ptr = “common part array; 
common part array = array {1..3] of common part rec; 


sub object rec = record 
num common parts : integer: 
common parts : common part ptr; 
sub _bsphere radius: real; {| RADIUS OF SUBOBJECTS BOUNDING SPHERE } 
sub bsphere x __: real; { CENTER OF BOUNDING SPHERE } 
sub bsphere y_ :: real; 


sub bsphere z_ : real; 
subobj type : Integer; { 0: SPHERE, 1: PLANAR-POLYGON } 
end; 


sub object ptr = “sub object array; 
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sub object array = array [1..3] of sub object _rec; 


object rec = record 
num sub objects : integer; 
sub objects : sub object ptr; 
obj bsphere radius: real; { RADIUS OF OBJECTS BOUNDING SPHERE } 
ob) bsphere x :real; { CENTER OF BOUNDING SPHERE } 
ob) bsphere y : real; 
ob) bsphere z : real; 
opcode : integer { CURRENTLY NOT USED } 


end; 


object ptr = “object array; 
object array = array [1..4] of object rec; 


light rec = record 


I ny Degqel -b»iereal; { INTENSITY OF THE LIGHT } 
light x, light y, light z:real; { LIGHT POSITION } 
dimension1, dimension? : real { NOT USED } 

end; 


light ptr = “light array; 
light array = array [1..3] of light rec; 


picture rec = record 


num objs : integer; 

objects : object ptr; 

num _ lights : integer; 

lights : hght_ ptr; 

global refraction index : real; 
no zero : real; 

ambient r : real; 

ambient g : real; 

ambient b : real; 


background color r : real; 
background color g: real; 
background color b: real; 


view position x_ : real; 

view position y_ : real; 

view position z_ : real; 

screen max X : integer; 

screen max y : integer; 
end: 


{ RKEKKKKKKEKKEEKAEKKEKEEEEEEKEKEREKKEKEKAKE EEE EKA EAE EEEKESE } 


raytype = (none, view, reflected, refracted); 
colortype = (red, green, blue); 


{ RAY DATA RECORD } 


ray ptr = “ray rec; 
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ray rec = record 


ray type : raytype; 
ray origin x :real; { ORIGIN OF RAY} 
ray origin y : real; 
ray origin z : real; 
ray vector x _:real; { DIRECTION OF RAY } 
ray vector y : real; 
ray vector z : real; 
ray stype -—-: raytype; { TYPE OF SOURCE RAY } 
intersection flag : boolean; 
obj idx ~: integer; { PATH DECSRIBING OBJECT INTERSECTED } 
subobj_idx : integer; 
cpart idx : integer; 
polygon idx : integer; 
intersection x :real; { INTERSECTION POINT } 
intersection y_ : real; 
intersection z_ : real; 
{ DISTANCE BETWEEN RAY’S ORIGIN AND INTERSECTION POINT } 
d : real; 


INTENSITY OF LIGHT COMING IN ALONG THE REFRACTED RAY GENERATED BY THIS 
RAY 


I tr, I tg, I tb : real; 


INTENSITY OF LIGHT COMING IN ALONG THE REFLECTED RAY GENERATED BY THIS 
RAY 


} 
I sr, 1 sg, | sb real, 
ray link : Tay ptr 
end; { ray ptr } 
EA 


{ KR KK KKK KKK KKK KAKA KEK KK KK AK KK KES VAR 2K oe oe eK ke OK KK KK KK KK KK KK KKK KKK KKK KKK K } 


var 


outfile : text; 
sysin : text; 
sysout : text; 


{ USED IN CONVERTING ALL RAYS/VECTORS TO UNIT VECTORS } 
x,y,z: real; 
NNILX, uNIty, uNnItz : real; 
jist : real: 


ray number : integer; 

{ USED AS SHORTHAND BECAUSE THE OBJECT PATHS GET LONG } 
cpart path : common part ptr; 
subobj path : sub object ptr; 


color: colortype; 
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{ USED TO IMPLEMENT THE STACK } 
ray top, 
ray current, 
ray next : ray ptr; 


intensity red, 
intensity zreen, 
intensity blue : real; 


pixel x, 
pixel y, 
pixel z : integer; 


ray generation number: integer; 


intersection point x, 
intersection point y, 
intersection point z_: real; 


old intersection point x, 
old intersection point y, 
old intersection point z : real; 


temp_integer|, 
temp integer2, 
temp integer3 : integer; 


reflected ray x, 
reflected ray y, 
reflected ray z : real; 


refracted ray x, 
refracted ray y, 
refracted ray z : real; 


surface normal x, 
surface normal y, 
surface normal z : real; 


{ USED TO INITIALIZE RAYS } 
initial ray type : raytype; 
initial ray origin x. initial ray origin y, initial rav origin 2 : real; 
imitial ray vector x. initial ray vector y, initial ray vector z : reali: 


— 


initial ray stype  : raytype: 
Initiai intersection flag : Dooiean; 
initial obj idx : integer; 
initial subobj idx : integer; 
initial cpart idx : Integer; 
initial polygon idx : integer; 


initial intersection x, 
initial intersection y, 
initial intersection z__ : real; 


ta 


initial d real; 
initial Bl mtr, en I tg, initial I tb: real; 
initial ‘Me _sr, initial I _ 8g, initial J _sb: real; 


{ HOLDS CURRENT RAY WHILE IN ACTUAL RAY TRACING LOOP } 
current ray type :raytype; 
current ray origin x, current ray origin y, current ray origin z : real; 
current ray vector x, current ray vector y, current ray vector 2: real; 


current ray stype : raytype; 
current intersection flag : boolean; 
current obj idx : integer; 
current _subobj | idx : integer; 
current cpart idx : integer; 
current polygon idx : integer; 


current intersection x, 

current intersection y, 

current intersection z_ : real; 

current d_ : real; 

current I _tr, current [ tg, current I tb: real; 
current Ts _sr, current fal _sg, current I _sb: real; 


reflected ray : boolean; 
refracted ray : boolean; 


source ray num : integer; 
source ray type: raytype; 


{ POINTER TO PICTURE RECORD } 


picture : picture rec; 


USED TO ESTABLISH LINKED LIST OF OBJ, SUBOBJ, LIGHTS, CPARTS, AND 
POLYGONS. 


} 
light _cntr : integer; 
light current : light ptr; 


ob) cntr : integer; 
ob) curr : object ptr; 


subob) cntr : integer; 
subobj curr : sub object ptr; 


cpart cntr : integer; 
cpart curr : common part ptr; 


poly cntr : integer; 
poly next, 
poly curr : polygon ptr; 


vertice cntr : integer; 
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{ USED WHEN READING IN DATA } 
poly loop cnt : integer; 
light loop cnt : integer; 
cpart loop cnt : integer; 
subobj loop cnt : integer; 
object loop cnt : integer; 
vertice loop cnt : integer; 


intersection X, intersection Y, intersection Z: real; 


{ USED TO IDENTIFY INTERSECTED OBJECTS } 
object idx, . 
subob) idx, 
cpart idx, 
polygon idx : integer; 


intersection flag : boolean; 


{ KRKKKRAKRKEKARAKEEKREKAKAKRAKHAKKEKKAKRHEKAKRKAKKRKAKKKAKRKAKAKAKAAKAKAKRAKRAKKAKKRKKAKAKEKXKE } 
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INTERSECTION PROCEDURES 


PESESSE SESE SE TESST ESE TES SSS ESS SE INTPRCS6 PAS KKK KKKKKKKKKKKKKKKKKKKKKKKEKKS 


* These are the intersection procedures used in the ray tracer program 
EKKKKKKKRKKKKKKKKKKKKKK KKK KKK KK KKKKKKKKKKRKKKKK KKK KKK KKKKKKKKKKKK KK KKKKKK KK 


} 


KEKKKKAKK KKK KKKKKKKKK KKK SPHERE INTERSECTION KKKKKKKKKKKKKKKKKHK KKK KKK KKK 

* CALLED FORM : CHECK FOR SUBOBJ-INTERSECTION 

* CALLS TO : NONE 

* DESC : CALCULATES THE INTERSECTION POINT BETWEEN A RAY/VECTOR AND A 
: SPHERE. 

* INPUT : The centerpoint of the sphere and it’s radius. The direction of 


bs the ray and a known point on the ray -- which would be it’s 

= origin. 

* OUTPUT : A flag indicating whether or not there was an intersection and 
y if there was the actual intersection point itself. 


KKKKKKKKKKKKKKKKKKK KK KAK EK KK KKK KK KK EK KEK KKKAKKKKKHKKKKAKRKKAKKKKAKKAKREKKKKKKKE 


} 


procedure sphere intersection (Px, Py, Pz: real; 
Vx, Vy, Vz: real; 
Cx, Cy, Cz: real; 
r : real; 
var o intersection flag : boolean; 
var Sphere x, Sphere y, Sphere z: real ); 


var 
dab cy bl, t2 : real; 
X1, Y1, Z1 : real; 
AZ, 12) 22 : real; 


distancel, distance2 : real; 

radical : real; 

diffx, diffy, diffz : real; 
begin 


{ INITIALIZE } 


Sphere x := 0; 
Sphere y := 0; 
Sphere z ak 
9 intersection flag := false; 


{ SET UP THE COMPONENTS OF THE QUADRATIC EQUATION } 
a := ( sqr(Vx) + sqr(Vy) + sqr(Vz)); 
b := ( (2 * Px * Vx) + 
(2 * Py * Vy) + 
(edge = 
(2 * Cx * Vx) - 
(2 * Cy * Vy) - 
(2 = eu Vz) ): 


74 


c := ( sqr(Cx) + sqr(Cy) + sqr(Cz) + 
sqr(Px) + sqr(Py) + sqr(Pz) - 
(25° Gx = Px) - (2 * Cy *#Py) - (2 * Cz * Pz) -*sqr(r)); 


radical := (sqrt ( sqr(b) - (4 * a * c))); 


{ START COMPUTATIONS ON QUADRATIC EQUATION } 
if (radical < 0) then begin 
writeln(’7 WARNING - imagionary number possible in SPHERE INTERSECTION’); 
o intersection flag := false 
end 


else begin 
if (radical = 0) then begin 


IF 0 THEN JUST ONE INTERSECTION POINT(THE LINE IS TANGENT TO THE SPHERE) 
} 
t1:=(-b:/ 2* a): { SOLVE FOR t } 
Sphere x := Px + (Vx * t1); { CALCULATE POINT USING t } 
Spitere y := Py — (Vy * t1); 
Sphere z := Pz + (Vz * t1); 
o intersection flag := true 
end 
else begin 


THERE WERE TWO INTERSECTION POINTS -- ONE ENTERANCE POINT AND ONE EXIT 
POINT. SOLVE FOR BOTH ts’ 
} 
tl : 


(-b + (sqrt ( sqr(b) - (4 * a *c)))) /2* a; 
2 "aa 


(-b - (sqrt ( sqr(b) - (4 *a*c)))) /2* a; 


X1 := Px + (Vx * t1); 
Y1 := Py + (Vy * t1); 
Z1 := Pz + (Vz * t1); 


CALCULATE DISTANCE FOR BOTH INTERSECTION POINTS FROM THE POINT OF ORIGIN 


distancel := (sqrt (sqr(X1 - Px) + sqr(Y1 - Py) + sqr(Z1 - Pz))); 
mezi= Px + (Vx * 2); 
heee:— Py + (Vy * t2); 
eee, (Vz * 2): 


dissance2 := (sqrt (sqr{X2 - Px)= sqr{Y2- Py) = sqr(Z2 - Pz))); 


{ COMPARE DISTANCES AND SELECT THE INTERSECTION POINT THAT IS CLOSER } 
if distancel < distance2 then begin 
Sphere x := X1; 
Sphere y := Y1; 
Sphere z := Z1; 
end 
else begin 
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Sphere x := X2; 

Sphere y := Y2; 

Sphere z := Z2; 
end; 


o intersection flag := true 


end: { ELSE} | 


diffx := Px - Sphere x; 
diffy := Py - Sphere y; 
diffz := Pz - Sphere 2; 


A CHECK TO INSURE THAT THE POINT SELECTED ISN’T THE RAYS POINT OF ORIGIN 


if ((diffx <= 0.0000) and (diffy <= 0.0000) and (diffz <= 0.0000)) then 


o intersection flag := false; 


end;{ ELSE } 
end; {* SPHERE INTERSECTIONS *} 


{PA} 


a CALCULATE PLANE EQUATION KKKKKKKAKAKKKKKKKKKKKK KK KK 
* CALLED FROM: FIND INTERSECTED POLYGON 

* CALLS TO : NONE 

* DESC : Calculates the constants of the equation of a plane when given 

; three point on the plane. 

* INPUT : Three vertices of a planar polygon. 

* OUTPUT : The A, B, C, and D constants for the equation of a plane. 


KAKKKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEK EK KKK KEKE KEKE KAKA KKK KEKE KEK KK 


} 


procedure calculate plane equation(X1, Y1, Z1 : real; 
X2AX 2,72 \: real; 
X3, Y3, Z3 : real; 
var A,B,C,D:: real ); 


begin 
A := Y1 * (Z2 - Z3) + Y2 * (Z3 - Z1) + Y3 * (Z1 - 22); 
B ss 7 (X2 eX) = 72 XX Zoe erie 
C := X1* (Y2 ~ VS ee) ee 


D :=-X1* ({Y2 * Z3) - (YB Z2)) omc) — (eee 
- (X3 * ((Y1 * Z2) - (Y2 * Z1))); 


end; { calculate plane equation } 
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KREKEKEKKEKKEEKEKKX KE FIND INTERSECTION POINT KREKKEKEKKKKEKKEKEKKKKE XS * 
* CALLED FROM: FIND INTERSECTED POLYGON, 
. CHECK FOR _SUBOBJ _INTERSECTION 
CHECK FOR INTERSECTION 
CALLS TO : NONE 
DESC : Calculates the intersection point between a ray/vector and 
a plane. 
INPUT : Ray direction and a known point on the ray, i.e., it’s 
point of origin. The constants (A,B,C,D) of the equation 
of a plane. 
OUTPUT : The intersection point. 


REKEKKEKKKEKKEKKK KECK EKEEKKKKEKKKKEKKEKE KEKEKKEKKKKEKKEKEKEKEKAKKKKKEKKEKEKKKK EX 


*% *% #=$8§ = 8 *% * 


procedure find intersection point(i A,i B,i C,i D : real; 
1 ray X,1 ray y,1 ray z : real; 
1 source x,1i source y,i source z: real; 
var o intersection point x, 
O intersection point y, 
O intersection point z : real); 


var 
t : real; 


begin 


{ SET UP FOR FINDING t FROM THE EQUATION FOR A LINE AND A PLANE } 
= (i_D- ((i_A *i_ source x) + (i_B*1i source y) + (i C * i source z))) / 
((i_A *i ray x) + 
(i B*i ray y)+ 
(i C *i ray z)); 


SUBSTITUTE t BACK INTO THE EQUATION FOR A LINE TO GET THE INTERSECTION 
POINT. 


} 


o intersection point x := (i ray x * t) +1 source x; 
o intersection point y := (i ray y * t) +1 source y; 
o intersection point z := (i ray z * t) +1 source 2; 


end; { find intersection point } 

{-PA} 

{ 

<M EREEKEEEEEKKLXx CALBO RRR MK RK SKE KE KK KE EKER KKK KEK KKK KKEK 
* CALLED FROM: FIND INTERSECTED POLYGON 

* CHECK FOR _SUBOBJ _INTERSECTION 

. CHECK FOR _INTERSECTION 

* CALLS TO : NONE 

* DESC: This calculates the constants A,B,C,D of the equation of a plane 
: given only a point on the plane and the surface normal of the 

plane. 

* INPUT : A point and the surface normal of the plane whose equation you 


* 


T7 


: are trying to figure out. 


* OUTPUT: The constants A,B,C,D of the equation of a plane. 


KAKA KKKKKKAKKAKKKKAAKKAKKAKKKAKKAKKKAKKKAKAKKAKKKAKRKKKKAAKRKKAKKKAKKK KKK KEK 


procedure caleq(bsphere x, bsphere y, bsphere z, 
vector x, vector y, vector z : real; 


var A, B, C, D: real); 


begin 
A := vector x; 
B := vector y; 
C := vector 2; 
D := (vector _x * bsphere x) + 


(vector y * bsphere_ y) + 
+ : 

’ 
(vector z * bsphere z) 


end; { caleq } 


\ hte eee POLYGON ORIENTATION KKKAKKEKKAKKKKKKKKKKKKKKKEKES 
* CALLED FORM: FIND INTERSECTED POLYGON 

* CALLS TO : NONE 

* DESC : This checks to see if the incoming ray will hit the front face of 

this polygon. It does this by comparing the angle between the 

surface normal of the polygon and the inverse of the incoming ray. 

If the angle is greater than 90 degrees then the ray is 

approaching the back of the polygon. 

* INPUT: The direction of the incoming ray. The surface normal of hte 

i object. 

* OUTPUT: A boolean value TRUE/FALSE depending on whether or not the 


_ polygon is facing the right direction. 
KAKAKKAKKKKAKKAKKKKKAKKKKAKKAKKKKAKKKKAKKKAKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKES 


} 


procedure polygon orientation (view vector x, view vector y, view vector z : real; 
surface normal x, surface normal y, surface normal 2: real; 
varo good orientation : boolean); 


* 
* 
* 
* 


var 
cosine theta : real; 
length view vector : real; 
length surface normal : real; 
dot product : real; 
begin 


TAKE THE DOT PRODUCT OF THE INVERSE OF THE VIEW VECTOR AND THE SURFACE 
NORMAL OF THE POLYGON IN QUESTION. 


} 
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dot product := (- view vector x) * surface normal x + 
(- view vector y) * surface normal y + 
(- view vector z) * surface normal 2; 


{ CALCULATE THE MAGNITUDE OF THE VECTORS } 
length view vector := sqrt(sqr(- view vector x) + 
sqr(- view vector y) — - 
sqr(- view vector z)); 


length surface normal := sqrt(sqr( surface normal x) + 
sqr( surface normal y) + 
sqr{ surface normal 2)); 


{ CALCULATE THE COSINE OF THE ANGLE BETWEEN THE RAYS } 


cosine theta := dot product / (length view vector * length surface normal); 


if (cosine theta > 0) then 

o good orientation := true 
else 

o good orientation := false; 


end; { polygon orientation } 


{ PA} 


RKEKKEKRAKEAKAEKKKKKKKKES FIND INTERSECTED POLYGON KAKKKKKKKKKKKKKKKKKKKEK KK ES 
* CALLED FROM : CHECK FOR SUBOBJ INTERSECTION 
* CALLS TO : POLYGON ORIENTATION 


: CALEQ 
: FIND INTERSECTION POINT 
‘ CALCULATE PLANE EQUATION 


* DESC : This determines if there is an intersection between a line/ray 

a and a polygon. If there is it calculates what it is. 

* INPUT: The direction of the shooting ray, a known point on that ray - 

: it’s origin, and the object path identifying the subobject to 

examine. 

* OUTPUT: A flag indicating whether or not a polygon was hit. If one was 
: hit the path identifying which one it was and the actual 


intersection point itself. 
KKK KEKE KAKA EE KKK EK EKA KKK KKK KKK KK AKERS 


} 


orocedure find intersected polygonii ray x,i ray y,i ray z: real: 


* 


* 


1 source X,1 source y,1 source z: real; 
1 obj idx, 
i subobj idx : integer; 
var o polygon intersection x, 
o polygon intersection y, 
o polygon intersection z : real; 
var o cpart idx, 
o polygon idx : integer; 
var o intersection flag : boolean); 
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type 


THIS IS SET UP TO HANDLE 6 SIDED POLYGONS. AS EACH SIDE OF THE POLYGON 
IS TESTED TO SEE WHETHER THE INTERSECTION POINT LIES INSIDE OR OUTSIDE OF 
IT THE CORRESPONDING ELEMENT IN THE ARRAY IS SET EITHER TRUE OR FALSE. AN 
ARRAY OF ALL TRUE MEANS THAT THE INTERSECTION POINT LIES WITHIN THE 
POLYGON, 
} 


intersection array = array(1..6] of boolean; 


var 
point outside polygon : boolean; 
intersections : intersection array; 


cpart cnt, 
polygon cnt, 
vertice cnt : integer; 


intersection found, 
good orientation : boolean; 


cpart path : common_part ptr; 
polygon path : polygon _ ptr; 


{ USED FOR ARBITRARY POINT TO DEFINE BOUNDING PLANE } 


anchor x, anchor y, anchor z: real; 


markerD, 
check point D: real; 


polygonX, polygonY, polygonZ : real; 


A, B, C, D: real; 
x1, yl, 21: real; 
x2, y2, 22: real; 
x3, y3, 23: real; 


diffx, diffy, diffz : real; 
loop cnt : integer; 
begin 


{ INITIALIZE } 
o cpart idx := 0; 
o polygon idx := 0; 
o polygon intersection x := 0.0; 
o polygon intersection y := 0.0; 
o polygon intersection z := 0.0; 
o intersection flag := false; 
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cpart cnt := 1; 

polygon cnt := 1; 

vertice cnt := 1; 

intersection found := false; 

good orientation := false; 

cpart path := picture.objects “{i obj idx|.sub objects *[i_subobj idx]. 
common parts; 

polygon path := cpart path “[cpart_cnt].polygons; 


THIS CHECKS EACH COMMON PART FOR AN INTERSECTION WITH ONE OF IT’S 
POLYGONS. 


} 


repeat 


{ THIS LOOP CHECKS EACH POLYGON OF A COMMON PART FOR INTERSECTION } 


repeat 
point outside polygon := false; 
polygon orientation(i ray x, 1 nay y,1 rey 2, 
polygon path “(polygon cnt|.surface normal x, 
polygon path “(polygon cnt].surface normal y, 
polygon path *[polygon_cnt| surface normal z, 
good orientation); , 


if good orientation then begin 


caleq(polygon path “(polygon cnt].vertice x(1], 
polygon path “(polygon cnt].vertice y(1], 
polygon path “|polygon cnt].vertice 2(1], 
polygon path “|polygon cnt].surface normal x, 
polygon path “|polygon cnt].surface “normal ys 
polygon path “[polygon cnt]. surface normal z, 
A,B,C,D): 


find intersection point(A, B, C, D, 
lyrey x; 
i ray_y, 
1 ray 2, 
i source x, i source y, i source 2, 
polygonX, Paiveonve polygonZ); 


CHECK TO MAKE SURE YOU ARE NOT CONSIDERING THE SOURCE POINT(ORIGIN) OF 
THE RAY. 


} 
diffx := 1 source x - polygonX; 
diffy := i source _y - polygonY; 
diffz := i source z - polygonZ; 


THIS SETS A FLAG IF YOU DO CONSIDER THE SAME POINT, THE CHECKS IN THE 
OTHER PROCEDURE SHOULD PREVENT THIS BUT JUST IN CASE. 
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if ((diffx <= 0.000) and (diffy <= 0.000) and (diffz <= 0.000)) 
then begin 


{ writeln(sysout,’set trip-wire’);} 
intersections/1] := false 

end 

else begin 


THIS LOOP CHECKS EACH EDGE OF THE POLYGON TO SEE IF THE INTERSECTION 
POINT LIES INSIDE OR OUTSIDE OF IT. 


} 
repeat 
{ THIS SELECTS THE FIRST VERTEX OF A POLYGON } 
x1:= polygon path [polygon cnt].vertice x[vertice cnt]; 
yl := polygon path “*|polygon cnt].vertice y|vertice cnt]; 
z1 := polygon path *[polygon cnt].vertice z[vertice cnt; 


if (vertice cnt = (polygon path “polygon _cnt/.num vertices - 


1)) 


WHEN YOU PICK THE NEXT TO LAST VERTEX YOU CAN SELECT THE NEXT 
CONSECUTIVE VERTEX TO ESTABLISH THE EDGE THROUGH WHICH YOU WANT THE 
BOUNDING PLANE TO PASS. YOU THEN MUST PICK ONE OF THE OTHER VERTICES (AND 
IT DOES NOT MAKE ANY DIFFERENCE WHICH ONE. I PICK THE FIRST ONE.) TO BE 
USED TO PUT INTO THE EQUATION OF THE PLANE THE RESULT OF WHICH IS COMPARED 
AGAINST THE RESULT THAT COMES FROM PLUGGING THE INTERSECTION POINT INTO 
THE EQUATION OF THE BOUNDING PLANE. 
j 
then begin 
x2 := polygon path “(polygon cnt]. 
vertice x[vertice cnt + 1]; 
y2 := polygon path “|polygon cnt]. 
vertice y|vertice cnt + 1]; 
22 := polygon path *[polygon cnt]. 
vertice z[vertice cnt + 1]; 


x3 := polygon path “(polygon cnt].vertice x(1]; 


y3 := polygon path “|polygon cnt|.vertice y(1]; 
z3 := polygon path “|polygon cnt|.vertice 2{1| 
end 
else begin 


f 
[F YOU DO®"NOT HAVE THE NEXT TO LAST EDGE WHEN Vig EEG Loe Naw 
CONSECUTIVE VERTICE TO ESTABLISH THE EDGE FOR THE BOUNDING PLANE AND 
THE ONE AFTER THAT TO PLUG INTO THE EQUATON OF THE PLANE. 
} 
x2 := polygon path [polygon cnt]. 
vertice x|vertice cnt + 1]; 
y2 := polygon path “|polygon cnt]. 
vertice y|vertice cnt + 1]; 
z2 := polygon path “(polygon cnt]. 
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vertice z|vertice cnt + 1]; 


x3 := polygon path “(polygon cnt]. 
vertice x|vertice cnt + 2]; 
y3 := polygon path “(polygon cnt]. 
vertice y{vertice cnt + 2]; 
23 := polygon path “'polygon cnt). 
vertice z(vertice cnt + 2| 
end; 


if (vertice cnt = polygon path “|polygon cnt|.num_vertices ) 


{ 
IF THE VERTICE SELECTED IS THE LAST ONE THEN JUST PICK THE FIRST VERTICE 
TO ESTABLISH YOUR BOUNDING EDGE AND THE SECOND VERTICE TO PLUG INTO THE 
EQUATION OF THE PLANE. 
} 
then begin 

x2 := polygon path *|polygon cnt].vertice x/1}; 

y2 := polygon path {polygon cnt].vertice y/1]; 

z2 := polygon path [polygon cnt|.vertice 2(1\; 


x3 := polygon path “[polygon cnt].vertice x(2]; 

y3 := polygon path “|polygon cnt].vertice y(2]; 

z3 := polygon path “polygon cnt|.vertice 2(2| 
end; 


{ 
ESTABLISH THE ARBITRARY POINT THROUGH WHICH THE PLANE WILL PASS. 
} 
anchor x := ((x2 + x1) / 2) + 10; 
anchor y := ((y2 + yl) / 2) + 10; 
anclior z := ((z2 + 21) / 2) + 10; 


calculate plane equation( x1, yl, z1, 
MEnOV 2, Voy 
anchor x, anchor y, anchor z, 


A, B, C, D); 


{ 

THE RESULT OF PLUGGING IN THE VERTICE OF THE POLYGON INTO THE EQUATION 
OF THE PLANE. 
} 


{ 
THE RESULT OF PLUGGING IN THE INTERSECTION POINT INTO THE EQUATION OF 


THE PLANE. 
} 


fener) <= (A ~ x3) + (Be ys) — (C473); 


PY 


check point D := (A * polygonX) + 
(B * polygonY) + 
(C * polygonZ); 
if ((markerD <= -D) and (check point D <= -D)) then 
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IF THE RESULTS HAVE THE SAME SIGN THEN THEY BOTH LIE ON THE SAME SIDE OF 
THE BOUNDING PLANE. HENCE THE INTERSECTION POINT LIES WITHIN THE POLYGON 
WITH RESPECT TO THAT EDGE. 
} 
intersections|vertice cnt] := true 
else 


IF THE RESULTS DON’T HAVE THE SAME SIGN THEN THEY LIE ON OPPOSITE SIDES 
OF THE BOUNDING PLANE. AT THIS POINT THE INTERSECTION POINT HAS BEEN 
PROVEN TO LIE OUTSIDE THE POLYGON. 


if ((markerD >= -D) and (check point D >= -D)) then 


intersections|vertice cnt] := true 
else 
intersections(vertice cnt] := false; 


vertice cnt — vertice cnt +1 


until (vertice cnt > polygon path “|polygon cnt|.num_ vertices); 
end 
end; 


CHECK THE POLYGON INTERSECTION ARRAY TO SEE IF THE INTERSECTION POINT 
FAILED THE INSIDE TEST FOR ANY OF THE EDGES. 
} 
for loop cnt := 1 to polygon path “[polygon_cnt].num_vertices do 
if not(intersections|loop cnt]) then 
point outside polygon := true; 


if point outside polygon then 
intersection found := false 
else 
intersection found := true; 


polygon cnt := polygon cnt + 1; 
until ((polygon cnt > cpart path *|cpart cnt].num_polygons) or 
(intersection found)); 


cpart cnt := cpart cnt + 1; 
until ((cpart_cnt > picture.objects “|i obj idx].sub objects “|i subobj_ idx]. 
num common parts) or (intersection found)); 


. SET UP THE OUTPUT POR THE PROGIPURE 
o cpart idx := (cpart cnt - 1); 
o polygon idx := (polygon cnt - 1); 
o polygon intersection x := polygonX; 
o polygon intersection y := polygonY; 
o polygon intersection z := polygonZ; 
o intersection flag := intersection found; 
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end; { find intersected polygon } 


{ PA} 


rESE SS ESSER SESE ES ES S | CHECK FOR SUBOBJ INTERSECTION HSHREKRKKEEHEEEEEEEKEH EE EXE XX 
* CALLED FROM: CHECK FOR INTERSETION 
* CALLS TO : CALEQ 


. FIND INTERSECTION POINT 
: SPHERE INTERSECTION POINT 
‘ FIND INTERSECTION POINT 


* DESC: Check to find out if the shooting ray intersects this 

¥ subobject’s bounding volume. 

* INPUT : The object whose bounding volume has been hit. 

; The direction of the shooting ray. 

The origin of the shooting ray. 

* OUTPUT : A flag indicating whether or not there has been an intersection. 
. The intersection point -- if there is one. 


The path to the intersected object. 
KREKKEKEKKKEKEKKEKRKEKKEEKKKE KEK KKEKKEKKEKEKEKKKKKEEKKKEKKKKKKEKKKKREEKEKEEEK KS 


} 


x 


x 


procedure check for subobj intersection (i object idx  : integer; 

jyray X,i rey y,1 ray z : real; 
1 S@Mirce x, 
1 source y, 
1 source z : real; 

var © intersection x, 
o intersection y, 
o intersection z : real; 

var o subobj idx, 
o cpart idx, 
o polygon idx : integer; 

var o intersection flag : boolean); 


var 
closest object : real; 


{ VECTOR BETWEEN THE RAYS ORIGIN AND INTERSECTION POINT } 
view polygon vector x, 
view polygon vector y, 
view polygon vector z : real: 


distance from intersection, 
distance from view position : real; 


subobj path: sub object _ptr; 
obj path : object_ptr; 


cpart cnt, 
polygon cnt, 
subobj cnt : integer; { USED TO GO INTO RESPECTIVE ARRAYS } 
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A, B, C, D: real: { CONSTANTS FOR EQUATION OF A PLANE } 


INTERSECTION POINT BETWEEN THE RAY AND PLANE THE BOUNDING CIRCLE IS 
INSCRIBED ON. 
} 

bplane intersection x, 

bplane intersection y, 

bplane intersection z : real; 


INTERSECTION POINT BETWEEN THE RAY AND A POLYGON. 
polygonX, polygonY, polygonZ : real; 


intersection flag : boolean; 


begin 


o subobj idx := 0; 

o cpart idx := 0; 

o polygon idx := 0; 

o intersection flag := false; 
o intersection x := 0.0; 

o intersection y := 0.0; 

o intersection z := 0.0; 


ESTABLISH A DEFAULT DISTANCE WITH WHICH THE ACTUAL DISTANCES WIL BE 
COMPARED. 
} | 

closest object := 10000.0, 

subobj path := picture.objects *[i_object idx].sub objects; 

ob) path := picture.objects; 


{ THIS LOOP CHECKS EACH SUBOBJECT OF AN OBJECT } 
for subobj cnt := 1 to obj path “[i_object_idx].num_sub_ objects do begin 


{ FIRST ESTABLISH THE PLANE ON WHICH TO DRAW THE BOUNDING CIRCLE } 
caleq(subobj path *!subobj cnt].sub bsphere x, 
subob} path “~'subod) cnt|.sub bsphere y, 
subobj path ~'subobj cnti.suo bsphere 2, 
1 Fay xX, 1 Fay y,i ray 2, 


ABC 2p): 


{ FIND THE INTERSECTION POINT ON THAT PLANE } 
find intersection point(A,B,C,D, 
1 ray X,1 ray y,i ray 2, 
1 source x,1 source y,1 source 2, 
bplane intersection x, 
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bplane intersection y, 
bplane intersection z); 


DETERMINE THE DISTANCE BETWEEN THE INTERSECTION POINT AND THE CENTER 
POINT OF THE BOUNDING SPHERE/CIRCLE. 
5 
distance from intersection := 
sqrt(sqr(subobj path “{subobj cnt].sub bsphere x - 
bplane intersection x) + 7 7 
sqr(subobj path “*|subobj cnt].sub bsphere y - 
bplane intersection y) + 
sqr(subobj path ~|subobj cnt].sub bsphere z - 
bplane intersection z) ); - 


{ IF THE DISTANCE IS LESS THAN OR EQUAL TO THEN YOU HAVE AN INTERSECTION } 
if (distance from _ intersection <= 
subobj path “|subobj cnt].sub bsphere radius) then begin 


if (subob) path“|subobj cnt|.subobj type = 0) then begin 


IF SUBOBJECT TYPE IS A SPHERE THEN YOU NEED TO USE THIS INTERSECTION 
PROCEDURE. 
} 
sphere intersection(i source x, i source y,i source z, 
1 ray X,1 ray _y,1 ray 2, 
subobj path “[subobj cnt].sub_ bsphere x, 
subobj path ~|subobj cnt].sub bsphere y, 
subob) path “|subobj cnt].sub_bsphere z, 
subobj path “{subobj cnt].sub_bsphere radius, 
intersection flag, 
polygonX, polygonY, polygonZ); 
cpart cnt := 1; 
polygon cnt := 0; 


end 
else begin 
{ IF IT ISN’T A SPHERE THEN USE THESE } 
find intersected polygon(i ray _x,1 ray _y,i ray 2, 
7 7 j_source_x, 1 source y,1 source Z, 
1 object idx, 
su bob} _cnt, 
polygonX. polygonY, polygondZ, 
cpart cnt, 
polygon cnt, 
intersection flag); 


end; { * ELSE * } 
if intersection flag then begin 


IF THERE HAS BEEN AN INTERSECTION THE ESTABLISH THE VECTOR BETWEEN THE 
ORIGIN OF THE RAY AND THE INTERSECTION POINT. 
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view polygon vector x := i source x - polygonX; 
view polygon vector _y :=1 source y - polygonY; 
view polygon vector z:=1 source z - polygonZ; 


CALCULATE THE MAGNITUDE (DISTANCE) OF THE RAY. 

distance from view position := sqrt(sqr(view polygon vector x) + 
sqr(view polygon vector y) + 
sqr(view_ polygon vector z)); 


» 


if distance from view position < closest object then begin 


COMPARE IT AGAINST THE OTHER DISTANCES AND SELECT THE CLOSEST ONE 
closest object := distance from _view_ position; 
o subobj idx := subob) cnt; 
o cpart idx :=cpart cnt; 
o polygon idx := polygon cnt; 
o intersection flag := intersection flag; 
o intersection x := polygonX; 
o intersection y := polygonY; 
o intersection z := polygond 

end 
end 
end 
end; {* FOR *} 


end; { check for object intersection } 


{PA} 


KKKKKKKKK KK KKK KK KE KKK EK CHECK FOR INTERSECTION KAKKKKK KK AERA KKK KKK KEES 
* CALLED FROM: MAIN AND CALCULATE INTENSITY 
* CALLS TO : CALEQ 
% FIND INTERSECTION POINT 
: CHECK FOR SUBOBJECT INTERSECTION 
* DESC : DETERMINES IF THERE IS AN INTERSECTION BETWEEN THE SHOOTING RAY 
- AND THE BOUNDING VOLUME OF AN OBJECT. 
* INPUT : Direction of the shooting ray. 
’ Origin of the shooting ray. 
Pointer into the object array. 
* OUTPUT : Flag indicating whether or not there was an intersection. 
7 If there was an intersection then then the intersection point. 
The path to the intersected polygon. 


KKK KK KK KKK KEKE EK KEK KEKE KK KKK KK KE KK KEKE EEK KEK KKK KEKE KEKE EERE KEKE KK EEE KE 


} 


* 


procedure check for intersection (i ray x,i ray y,i ray z: real; 
1 source x,1 source y,1 source z: real; 
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i picture objects : object ptr; 
var oO intersection x, 

o intersection y, 

o intersection z : real; 
var o object idx, 

o subobj idx. 

o cpart idx, 

o poly idx : integer; 
var o intersection flag : boolean); 


var 
closest object : real; 


object cnt, 

object idx, 

subobj cnt, 

cpart cnt, 

polygon cnt : integer; 


distance from viewposition, 
distance from intersection : real; 


INTERSECTION POINT BETWEEN SHOOTING RAY AND THE PLANE THE BOUNDING CIRCLE 
IS INSCRIBED ON. 


} 
bplane intersection x, 
bplane intersection y, 
bplane intersection z : real; 
view bplane vector x, 
view bplane vector y, 
view bplane vector z : real; 
INTERSECTION POINT BETWEEN THE SHOOTING RAY AND A POLYGON. 
polygonX, polygonY, polygonZ : real; 
CONSTANTS FOR THE EQUATION OF A PLANE. 
} 
#, B, C, D eal: 
found intersection : boolean; 
distance : real; 


begin 


o object idx := 0; 
o subobj idx := 0; 
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o cpart idx := 0; 

o poly idx := 0; 

o intersection x := 0.0; 

o intersection y := 0.0; 

oO intersection 2 := 0.0; 

o intersection flag := false; 
closest object := 10000.0; 


found intersection := false; 


object cnt := 1; 


LOOP TO CHECK EACH OBJECT IN THE OBJECT ARRAY. 
} 


repeat 


if (((picture.objects “[object_cnt].obj bsphere z >i source z) and 
(i ray z>0)) or 
((picture.objects “jobject cnt|.obj) bsphere z <1 source z) and 
(i ray z<0))) then begin 


ESTABLISH PLANE ON WHICH TO DRAW BOUNDING CIRCLE. 
} 
caleq(picture.objects “|object cnt].obj bsphere x, 
picture.objects “|object_cnt].obj bsphere y, 
picture.objects “[object_cnt].obj bsphere z, 
1 ray x,1 ray y,1 ray 2, 


A, B,C, 


FIND INTERSECTION POINT BETWEEN THAT PLANE AND THE SHOOTING RAY. 
} 
find intersection point(A,B,C,D, 

1 ray X,1 ray y,i ray 2, 

i source x, 1 source y,1 source 2, 

bplane intersection x, 

bplane intersection y, 

bplane intersection z); 


DETERMINE DISTANCE BETWEEN THE CENTER OF THE CIRCLE AND THE 
INTERSECTION POINT. 
j 
distance from intersection := 
sqrt(sqr(picture.objects “[object_cnt|.obj bsphere x - 
bplane intersection x) + 
sqr(picture.objects “|object_cnt].obj bsphere y - 
bplane intersection y) + 
sqr(picture.objects “[object cnt].obj bsphere z - 
bplane intersection 2) }; 
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IF INTERSECTION POINT LIES WITHIN CIRCLE THEN START CHECKING THE 
SUBOBJECTS THAT MAKE UP THE OBJECT. 


if distance from intersection <= 
picture.objects “lobject cnt|.obj bsphere radius then begin 


check for subobj intersection (object cnt, 
1 ray X,1 ray y,i ray 2, 
1 source X,1 source y,1 source 2, 
polygonX, polygonY, polygonZ, 
subob) cnt, 
cpart cnt, 
polygon cnt, 
found intersection); 


DETERMINE THE DISTANCE BETWEEN THE ORIGIN AND INTERSECTION POINT OF 
OF THE RAY 


if found intersection then begin 
view bplane vector x :=1 source x - bplane intersection x; 
view bplane vector y :=1 source y - bplane intersection y; 
view bplane vector z:=1 source z- bplane intersection 2; 


distance from _viewposition := sqrt(sqr(view_bplane vector x) + 
sqr(view bplane vector y) + 
sqr(view _bplane vector z)); 


SELECT ONE CLOSEST TO RAY’S ORIGIN. 
j 

if distance from viewposition < closest object then begin 
closest _object := = distance from _viewposition; 
o object idx := object | cnt; 
Oo _subobj | ie j= subobj ent; 
o cpart idx := cpart cnt; 
o poly idx := polygon cnt; 
o intersection x := polygonX; 
o intersection y := polygonY; 
o intersection z := polygonZ; 
o intersection flag := found intersection 

end 

end 
end; 


end; 
object cnt := object cnt + 1; 


until (object cnt > picture.num_objs); 


end; { check for intersection } 
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INTENSITY PROCEDURES 


KKK KK KK RK KKK KKK KK KKK KKK KEKE KE STACK EMPTY KKK AK KKK KKK KK KKK KKK KKK KK EK 


* CALLED FROM : MAIN and POP 

“CALLS@O :NONE 

* DESC : CHECKS TO SEE IF STACK EMPTY 
* INPUT : POINTER TO TOP OF STACK 

* OUTPUT : BOOLEAN VALUE - TRUE/FALSE 


OK KKK RE KKK KR KKK KKK KKK KKK RK KKK KKK KEK KKK KK KKK KKK KK KKK KKK KK EKA KEK KKK KEK EKE KES 


} 


function stack empty (input ray top: ray ptr) :boolean; 
begin 
STACK EMPTY := input ray top = nil; 


end; { STACK EMPTY } 


KKK KK KEK KK KA KKAEKAEKAEKKEKKEK KE STACK EXCEEDED KKKKKKKKKKKKKKKKAKKKK KK KKE KS 


* CALLED FROM : MAIN 
* CALLS TO : NONE 
* DESC  : CHECKS TO SEE IF STACK IS FULL 
*INPUT :CURRENT SIZE OF STACK and MAXIMUM SIZE OF STACK 
* OUTPUT : BOOLEAN VALUE TRUE/FALSE 
kK Ro kK KR EK KK KK KOK KR KKK KKK KKK KK KK KR KK KK KK KK KKK KKK KK KK KKK KKK KKK RK KK KK KKK 
} function stack exceeded{input tl, 
input t2: integer) : boolean; 


begin 
if input tl > input t2 then 
stack exceeded := true 
else 
stack exceeded := false; 


end STACK EXGEED ie 


{PA} 


cae CALCULATE REFRACTED RAY KHKKKKKKKKKAKKKKEKKKKKKKKKK SE 


* CALLED FROM : MAIN 
* CALLS TO : NONE 
* DESC _ : This calculates the direction of a refracted ray. 
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* INPUT: Direction of source ray, surface normal of intersected object, 
s the objects index of refraction, the global index of refraction. 

* OUTPUT : A flag set to true or false, depending on weither or not a 

i refracted ray was created. If one was then it’s direction is 

given. 


‘ESSE SESS SESE ESSE EEE SE SESS SSS SESE SSS SSS SSS SSL SSS SESS SSS SSE SESS SSS SS SSS ES SS FS 


} 


procedure calculate refracted ray (i ray vector x, 

1 ray vector y, 
1 ray vector z : real; 
1 ob) surface normal x, 
1 obj surface normal y, 
1 ob) surface normal 2: real; 
1 obj ridx : real; 
1 global ridx =: real; 

var o refracted ray x, 
o refracted ray y, 
o refracted ray 2: real; 

var o refracted ray flag : boolean ); 


x 


var 
test Kf, 
testKf2, 
Kf : real; 
Kn : real; 
abs dot product : real; 


length of ray : real; 
VIExi vl y;vi.z : real; 


begin 


THIS PRODUCES THE ABSOLUTE VALUE OF THE DOT PRODUCT OF THE INCOMING RAY 
AND THE SURFACE NORMAL OF THE INTERSECTED SURFACE. 
} 
abs dot product := Abs((i obj surface normal x *i ray vector x) + 
(i obj surface normal y * 1 ray vector y) + 
(i obj surface normal z * i ray vector z)); 


if abs dot product = 0 then begin 


THIS [S JUST A PRECAUTION. 
o retracted ray dag := false; 
o refracted ray x := 0; 
o refracted ray y := 0; 
o refracted ray z :=0 
end 
else begin 


THIS PRODUCES THE UNIT NORMAL VECTOR IN THE DIRECTION OF THE INCOMING 
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ee 


} 
vl x :=1_ray vector x / abs dot product; 
vl_y :=i ray vector y / abs dot _product; 
vl z :=i_ray_ vector z / abs dot product; 
length of ray := (sqrt(sqr(vl x) ~ sqr(vi_y) + sqr(vl_z))); 


{ THIS GIVES THE RATIO OF THE REFRACTIVE INDICES } 
Kn := i obj ridx / i global ridx; 


THE CALCULATION OF THE FRESNEL COEFFICIENT IS DIVIDED UP THIS WAY 
INORDER TO CHECK LATER FOR AN IMAGINARY DENOMINATOR WHICH INDICATES 
TOTAL INTERNAL REFLECTION. 

} 
testKf := (sqr(Kn) * sqr(length of ray)); 
testKf2 := sqr(sqrt(sqr(i_ obj surface normal x + v1_x) + 
sqr(i_obj surface _normal ile vi oy lea 
sqr(i_ obj : “surface “normal 12 3. Viz) 


if (testKf - testKf2 <= 0) then begin 
{ IMAGINARY DENOMINATOR - TOTAL INTERNAL REFLECTION IS OCCURING } 
o refracted ray flag := false; 
Oo _refracted Tay x := 0; 
Oo _refracted _Tay_ y:=0,; 
o refracted ray z: —O 
end 
else begin 
o refracted ray flag := true; 
Kf := 1 / sqrt(testKf - testKf2); 
o refracted ray x := (Kf * (i_obj surface normal x + v1_x)) - 
1 obj surface normal x; 
o refracted ray y := (Kf * (i obj surface normal y + v1 y)) - 
1 ob) surface normal y; c 
o refracted ray z:= (Kf * (i obj surface normal z+ v1 z)) - 
1 obj surface normal z 
end 


end; 
end: { calculate refracted ray } 


ea 

i 

We KR OK ok ke kk kk kk ok ok ok ok ok ok ok ok CALCULATE REFLECTED RAY 3 kk kk ok ok kok ok ok ok ke ok ke ok ok kk ok ok ok 
* CALLED FROM: MAIN and CALCULATE INTENSITY 

* GALLS TO +: NONE 

*DESC  : CALCULATES THE DIRECTION OF A REFLECTED RAY 

*INPUT : Direction of the source ray and surface normal of intersected 

- object. 

* OUTPUT : Flag indicating existance of reflected ray. If one exists 
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< then it’s direction is given. 


KEK KKK KKK EK ERK KKK KK KEKE EK KKK KEK KKK KKK KKK KK EK EK KEKE KK KEKE KKK KEK EKAKKEKKK EEK 


procedure calculate reflected ray (i ray vector x, 
1 ray vector y, 
1 ray vector z : real; 
1 obj “surface _normal x, 
1 obj. “surface _normal | oN 
1 obj : “surface “normal Tae teal; 
var o reflected ray x, 
o reflected ray _y, 
o reflected Mray. 2’: real: 
var o reflected _ray flag : boolean ); 


var 
abs dot product : real; 
length | of ray : real; 
Tie v1, v1 2 : real 
Ri * RL ae Rl wae neal; 


begin 


ABSOLUTE VALUE OF THE DOT PRODUCT OF THE OBJECTS SURFACE NORMAL AND 
OF THE INCOING LIGHT RAY. 
} 
abs dot product := Abs((i_ obj surface normal x *i ray vector x) + 
(i_obj surface normal y * 1 ray_ vector -y) + 
(i_obj surface normal z *i ray vector z)); 


vl_x :=i ray vector x / abs dot product; 
vl_y :=i ray vector y / abs dot product; 
vl _z :=1i ray vector z/ abs dot product; 


{ CALCULATE UNIT NORMAL VECTOR IN THE DIRECTION OF THE INCOMING RAY. } 
Ri x :=vil_x + (2*i obj surface normal x); 
Rly := vl_y + (2 *1i obj surface normal y); 
Ri z == vl _z.+ (2 * i obj surface normal z); 


length of ray := (sqrt(sqr(R1_ x) + sqr(R1_y) + sqr(R1_z))); 
{ CALCULATE REFLECTED RAY } 
o reflected ray x:= R1 x, length of ray; 
o reflected ray y:=R1 y_ length of ray; 


Oo _reflected vay 2:=R1 _z / lengtn of ray; 
Oo “reflected _Tay Wikg := true; 


end; { calculate reflected ray } 


{ PA} 
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KKKKKAKKKKKAKKKKKKKKKKKKKS CALOCUPTATE INTENSITY KHKKKKKKHKKKKKKKKKKKKKKS 


* CALLED FROM : [2 and I6 
* CALLS TO 
* DESC 

* INPUT 


= 


; OUTEW 


- set in the appropriate source ray. 


: CHECK FOR_INTERSECTION and CALCULATE REFLECTED RAY 
: Calculates the intensity at any given point. 
: Color component being calculated, current ray data, 
pointer to object data, and pointer to light data. 
: Intensity at a given point, either to be displayed or 


KKKKAKKAKKKAKKKKKKKK KK KKKAKKKKKKKKKKKK KKK KKKKAKKKHKKKAKAKKKAKKKKKKKKKKKKE KS 


procedure CALCULATE INTENSITY ( input color 


input ray d : Teal; 
input ambient : real; 
input ray I t : real; 
input ray I s : real; 


input ray vector x, 
input ray vector y, 


input ray vector z : real; 


: colorty pe; 


input | number _of lights : integer; 
input “obj _ K a : Teal; 

input obj __ ie s : real; 

input obj | “KS t : real: 

input obj _ K _d : real; 

input obj phong exp : integer; 


input intersection x : real; 
input intersection y : real; 
input intersection 2 : real; 
input _obj surface “normal _X: 
input obj : “surface normal - oye 
input obj | surface “normal 2 
input picture object : object 
input light top : light ptr; 
var io intensity : real }; 


var 
ans : 
I.d 
I | 
j, 1: integer; 
distance : real; 
sight ray x, 
sight ray 3 


char; 
: real; 
sreal; 


sight ray z . real: 
unit “sight X, 
unit sight y, 
unit sight z 
light ray x, 
light ray y, 
light ray z 

unit light | x 
umit light y, 


: real; 


: real: 
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real; 
real; 


: real; 


_Ptr; 


unit light z : real; 

unit reflected mx, 

unit reflected — S¢ 

unit reflected z : real: 
intersection x, 
intersection y, 

intersection 2 : real; 
reflected _light eray xX, 
reflected light ray y, 
reflected light ray z : real; 
source X, 

source y, 

source z : real; 

check x, 

check y, 

check z : real; 

obj idx, subob) idx, cpart idx, polygon idx : integer; 
intersection flag : boolean; - 
ob) light distance : real; 
reflected ray : boolean; 
intersected obj Kt: real; 


begin 


{ THIS ELIMINATES THE SAME INTERSECTION POINT FROM BEING SELECTED AGAIN.} 


if input ray d > 0.1 then begin 
I 1 := 0.0; 
io intensity := 0.0; 
reflected ray := false; 
intersection flag := false; 


{ THIS SETS UP THE SIGHT RAY } 


sight ray x :=-input ray vector x; 
sight ray _y :=-input ray vector y; 
sight ray z:=-input ray vector 2; 


distance := sqrt(sqr(sight _ray_x) + 
sqr(sight ray y) + 
sqr(sight ray z)); 


unit sight x := sight ray x / distance; 
unit sight v := sight ray y / distance: 
unit sight z:= signt ray z distance; 


fori := 1 to (input number of lights + 1) do begin 
{ THIS GENERATES THE SHADOW FEELERS } 


light ray x := picture.lights “[i].light x - input intersection x; 

light ray y: picture.lights “[i].light_y - input intersection _y; 

light ray 2 := picture.lights *[i].light_z- input intersection 2; 

distance := (sqrt(sqr(light ray x) + 7 
sqr(light ray y) + 
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sqr(light ray z))) 


{ CONVERTS IT TO A UNIT VECTOR } 
unit light x := light ray_x / distance; 
unit light y := light ray y / distance; 
unit light z := light ray z / distance; 


source x := picture.lights “|i].light x; 
source y := picture.lights “[i].light_y; 
source z := picture.lights “[i].light _z; 


CHECK TO SEE IF ANY OF THE SHADOW FEELERS INTERSECT ANYTHING. 
} 
check for intersection(-unit light x, -unit light y, -unit light z, 

- source X, 
source y, 
source 2, 
Input picture object, 
intersection x, 
intersection y, 
intersection 2Z, 
obj idx, subob) idx, cpart idx, polygon idx, 
intersection flag); 


CHECK TO INSURE THAT THE SAME POINT IS NOT CONSIDERED AGAIN, WHICH CAN 
HAPPEN. 


} 


check x := intersection x - input intersection x; 
check y := intersection ‘y - input intersection y; 
check z := intersection z- input intersection Zz; 


PULL THE PROPER CHARACTERISTICS OF THE OBJECT OUT TO DEAL WITH THE 
APPROPRIATE COMPONENT OF LIGHT THAT IS BEING CURRENTLY DEALT WITH 
} 
if color = red then 
I ]1:= picture.lights “(i].I_r; 
intersected obj Kt := picture. 
objects “[obj_ idx]. 
sub objects ~'subob} idxl. 
common parts icpart 1dxj. 
Nictr, - - 
if color = green then 
I 1 := picture.lights “[i].I_g; 
intersected obj Kt := picture. 
objects *{obj idx]. 
sub objects “|subobj idx]. 
common parts *[cpart idx]. 
K tg; 
if color = blue then 
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I_1:= picture.lights “li].I_b; 
intersected obj Kt := picture. 
objects *[obj _idx]. 
sub objects *{subob} _idx]. 
common parts “[cpart idx]. 


K_ tb; 


obj light distance := (sqrt(sqr(intersection x - 
picture.lights “[i].light x) + 
sqr(intersection y - 
picture.lights “[i].light y) + 
sqr(intersection z - - 
picture.lights “|i].light_z))); 


IF THERE HAS BEEN AN INTERSETION AND THE POINT BEING CONSIDERED IS NOT 
THE ORIGIN OF THE RAY THEN CHECK TO SEE IF THE OBJECT IS OPAQUE. IF IT IS 
THEN RETURN TO BEGINNING OF LOOP AND CHECK NEXT SHADOW FEELER. IF IT IS 
NOT OPAQUE THEN CALCULATE INTENSITY AT THAT POINT AND CONTINUE FOLLOWING 
THE RAY TO SEE IF IT INTERSECTS ANYTHING ELSE. CONTINUE THIS LOOP 
EITHER UNTIL NO MORE OBJECTS ARE LEFT OR UNTIL AN OPAQUE ONE IS 
INTERSECTED. 
| 
THIS PART OF THE CODE HAS NEVER BEEN TESTED. I JUST CODED IT AS I 
THOUGHT IT SHOULD BE FROM THE ALGORITHM IN ROGERS BOOK PP. 377. 
QUITE FRANKLY I STILL DON’T FULLY UNDERSTAND WHAT IS SUPPOSE TO TAKE PLACE 
HERE. 
} 
if intersection flag then begin 
while ((intersection flag) and 
((check x > O) or 
(check y >) or 
(check z > 0))) do begin 


intersection flag := false; 


if not( intersected obj Kt = 0) then begin 
if (input color = red) then begin 


I 1:= picture.lights “[i].I r * 
- intersected obj Kt: 
end; 
if (inpuc color = green) then begin 
I 1:= picture.lights “[i].I g * 
intersected obj Kt; 
end; 
if (input color = blue) then begin 


I 1:= picture.lights *[i].I b * 
intersected obj Kt; 
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end: 


calculate reflected ray(input ray vector x, 
input ray vector y, 
input ray vector z, 
input obj surface normal x, 
input obj surface normai y, 
input obj surface normal 2, 
reflected ray x, 
reflected ray y, 
reflected ray z, 
reflected _ray); 


distance := sqrt(sqr(reflected ray x) + 
sqr(reflected ray _y) + 
sqr(reflected ray z)); 


unit reflected x := reflected ray x / distance; 
unit reflected | _y := reflected ray y / distance; 
unit _reflected _z:= reflected ray z / distance; 


10 intensity := io intensity + 
(1.1 * input obj K d) * 
((input _obj : surface _normal x * unit light x) + 
(input obj | “surface “normal _y * unit light y) + 
(input obj surface normal z * unit light z))) + 
(ie input | _obj K _s) 7 
((unit sight x * unit reflected x) + 
(unit sight y * unit reflected y) + 
(unit sight z * unit reflected z))); 


source xX := intersection x; 
source y := intersection y; 
source Z := intersection 2; 


check for intersection(-unit light x, 
-unit light y, 
-unit hight z, 
source x, 
source y, 
source 2, 
input picture object, 
intersection X, 
intersection y, 
intersection 2, 
obj idx, subobj idx, 
cpart idx, polygon idx, 
intersection flag); 


intersection x - input intersection x; 


check x: 
check y := intersection y - input intersection _y; 
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check z := intersection z- input intersection 2; 


end {if} 
end {while} 
end {if} 
else begin 
calculate retlected ray({input ray vector x, 
input ray vector y, 
input _Tay vector 2, 
input obj ; “surface normal me 
input obj surface normal y, 
input ob) surface normal z, 
reflected ray x, 
reflected ray y, 
reflected ray z, 
reflected ray); 


distance := sqrt(sqr(reflected ray x} + 
sqr(reflected ray y) + 
sqr(retiected ray 2z)); 


unit reflected x := reflected ray x / distance; 
unit “reflected _ _y :=reflected ray y / distance; 
unit reflected 2 := reflected ray z / distance; 


lo intensity := lo intensity + 
((I1_1 * input obj K_d) * 
((input obj. surface _normal x * unit light x) + 
(input obj : “surface “normal _y * unit light y) + 
(input obj | surface normal 2 * unit light z))) - 
((I1_1 * input obj K 's) * 
((unit_sight_x * unit reflected x) + 
(unit sight y * unit reflected _y) + 
(unit sight z * unit reflected z))); 
end { ** ELSE **} 
end; { ** FOR ** } 


{ 
THIS IS THE STUB TO JUST HAVE EVERY OBJECT ILLUMINATED BY AMBIENT LIGHT 


io intensity := input obj K a * input ambient; 


{ RKAKKEKEKKKAKEKKKEKEKEKEKEEEKEKKKEKEKEKEKKEEKEEKEKEKEKKEKEK EK KK \ 


THIS IS WHERE THE FINAL INTENSITY IS CALCULATED 
S 
io intensity := ((input ob) K a * input ambient) + 
lo intensity + 
((input _obj K s* input ray I s) + 
(input obj K t * input nay J Cie 
{/input_ray dor /2 or /1} 
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THIS IS JUST TO KEEP ALL VALUES WITHIN A RANGE WHERE THEY CAN BE DISPLAYED 
} 
if io intensity > 1.00 then 
io intensity := 1.00; 
if io intensity < 0.00 then 
lo intensity := 0.00; 


{ ee 2 Eo 2 EO EK KK OOK Ek \ 


end 


else begin 


IF THE INPUT RAY DISTANCE IS LESS THEN ONE THAN YOU ARE CONSIDERING THE 
SAME POINT AND HENCE THE INTENSITY THERE SHOULD BE 0. 
} 

io intensity := 0.0 
end; 


end; { CALCULATE INTENSITY } 


{PA} 


ok kK eK OK KK KK OK KKK KEK KK KK KKK DISPLAY PIXEL OE eR eK KK kK eK ok OK KK KOK KOK KK KE 


* CALLED FROM : MAIN 

* CALLS TO : NONE 

* DESC: WRITES OUTPUT TO FILE 

* INPUT : THE VALUES FOR THE RED, GREEN, AND BLUE COMPONENTS OF LIGHT. 
* OUTPUT : NONE 


2 ok 2K oo ok oR ee ko 2 oe 2 Kk oe kK OK ke oO ok ok 2K 2K ok 2 OK ok ok KK OK KK kK KK KK KK KK KK KE KK KK KKK 


} 

procedure DISPLAY PIXEL ( input intensity red, 
input intensity green, 
input intensity blue : real; 
input pixel x, 
input pixel y, 
input pixel z : real }; 


begin 


write (outfile,input intensity red:3:2); 
write (outfile,’ ’,input intensity green:3:2); 
writeln (outfile,’ ’,input intensity blue:3:2); 


end; { DISPLAY PIXEL } 
clea 


ok 2K kk ok OK kK OK OK KKK KOK KOK KK KK OK KK KK POP eK KK ok OK KK ok kK OK KO KK KK OK KKK KK KKK KKK KKK KKK 


* CALLED FROM : MAIN, I2, I6 

* CALLS TO : NONE 

* DESC : Removes a ray from the top of the stack. 

* INPUT : Pointer to the current top top of stack. 

* OUTPUT : The ray just popped from the stack and a pointer to the new top 
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* 


top of stack. 


ORE KK KK OR OK KK OK OR OR KR OK OK KKK OK OK EK OK OK KK OK OK KK KK OK OK OK KK KK KOK OK OK OK OK OK KOK KK OK KKK KKK 


} 


procedure POP (var output ray type : raytype; 
var output ray origin x : real; 
var output ray origin y : real: 
var output ray origin z : real; 
var output fay vector x : real; 
var Output ray vector y : real; 
var output ray vector z : real; 
var output ray stype > raytype; 
var output intersection flag : boolean; 
var output obj idx oy integer; 
var output subobj idx : Integer; 
var output cpart idx : Integer; 
var output polygon idx : Integer; 
var output intersection x_ : real; 
var output intersection y_ : real; 
var output intersection 2 : real: 
var output d : real; 


begin 


var output I tr, output I tg, output I tb : real; 
var output I sr, output I sg, output I sb : real; 
var io top : ray ptr ); 


if (stack empty (io top)) then begin 
writeln(’STACK UNDERFLOW ERROR’) 


end { if } 

else begin 
output ray type := 10 top”.ray type; 
output ray origin x := 10 top”.ray origin x; 
output ray origin y = 10 top”.ray origin y; 
output ray origin z [= 10 top”.ray origin 2; 
output ray vector x = 1o_top “Pay vector x; 
output ray vector y := 10 top”.ray vector y; 
output ray vector z = 1o_top “ray vector 2; 
output ray stype := 10 top”.ray_stype; 
output intersection flag := io top ~.intersection flag; 
output obj idx := 10 top”.obj idx; 
output subobj idx := 10 top’.subobj idx; 
output cpart idx = A0 “top “cpart idx: 
outpus polygon idx = io top ’.polygon dx; 
output intersection x ‘= 10 tOp intersection x: 
output intersection y := 10 top”.intersection y; 
output intersection z := 10 top”.intersection 2; 
output d := 10 top .d; 
output I tr :=1o top’.I tr; 
output I| tg :=1o top”.] tg; 
output I tb :=1o top”.I tb; 
output I sr := lo top”.I sr; 
output I sg := 10 top”.] sg; 
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output [ sb := 10 top~.I sb; 
ray next := 10 top”.ray link; 
{ remove old pointer } 
_ dispose(io top); 
{ set top of stack pointer to new top of stack } 
lo top := ray next; 
end; { else } 


end one 


{PA} 


KAKKKKKKKKKKKKKKKK KKK KKKAKK KKK PUSH KK KKK KKK KKK KKK KKK KKEKKKKKKKKKKES 


* CALLED FROM: MAIN, 12, 16 

* CALLS TO : NONE 

* DESC : Places a ray on the top of the stack. 

* INPUT: The current top of the stack, and the data for a new ray 
* OUTPUT: The pointer to th new top of stack. 


2 KK KK KKK KK KK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKKKKKKKKKKKKKKKKKK 


} 


procedure PUSH ( input ray type : raytype; 
input ray origin x : real; 
input ray origin y : real; 
input ray origin z : real; 
input ray vector x : real; 
input ray vector y : real; 
input ray vector z : real; 
input ray stype : raytype; 
input intersection flag : boolean; 
input obj idx : integer; 
input subobj idx : Integer; 
input cpart idx : integer; 
input polygon idx : integer; 
input intersection x_ : real; 
input intersection y_ : real; 
input intersection z_ : real; 
input d : real; 


input I tr, input | tg, input I tb: real; 
input I sr, input I| sg, input I sb: real; 
var io top : ray ptr ); 


begin 


mew(ray current}; 


ray current .ray type = input ray type; 

ray current .ray origin x ‘= input ray origin x; 
ray current .ray origin y ‘= input ray origin y; 
ray current .ray origin 2 := Input ray origin 2; 
ray current “.ray vector x = input ray vector x; 
ray current "ray vector y = input ray vector y; 
ray current “.ray vector z = input ray vector 2; 
ray current .ray stype := Input ray stype; 


ray current’ .intersection flag := input intersection flag; 
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ray current ~.obj idx := input obj idx; 


ray current ~.subob) idx := input subobj idx; 
ray current ”.cpart idx := input cpart idx; 
ray current .polygon idx := input polygon idx; 
ray current .intersection x := input intersection ne 
ray current .intersection y := input intersection y; 
Tay current .intersection Z := input intersection 2; 
ray current °.d := input d; 
ray current”. tr := input I tr; 
ray current’.] tg = input I tg; 
ray current”.J tb := input I tb; 
ray current e .aer := input I sr; 
ray current’. sg := input I sg; 
ray current ~.I sb := input I sb; 
ray current .ray link [= 10 top; 
1o_top := Tay current; 

end; { PUSH } 
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MAIN 


PROGRAM RAYTRACER, 


KKK KK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKKRKKKKR KKK AKKKKAKAKKKK KES 


ERK MM KKK KK KKM KK KK KM KKK KR KK KK KE KKK KKK KKK KKK KKK KK KKK KKK KK KKK KK KK 


“SrROG > RAY.PAS 

** AUTHOR : Paul G. Smith 

** DATE : 11 May 1987 

** DESC : A ray tracing prototype with a global illumination model 


se integrated into it. 


** INPUT : A sequential scene file under PICTURE.PAS. 
** OUTPUT: A bitmap file containing the red, green, and blue color 
oi color components. Their values range from 0-1 and need to 


a be converted for display on an RGB color monitor. 
KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKKKKKKKKKKKKKKKKKKKSE 


KKKKKKKKKAK KKK KAK KKK KKK KKKA KK KKK KKK KKK KKK KKK KKK ARKKKKK KKK AKRKKKAKKKKKKKEKE 


} 
{* INCLUDE FILES *} 


{$I declare6.pas} - { DECLARATION SECTION } 

{.PA} 

{$I intpres6.pas} - { INTERSECTION PROCEDURES } 

eA} 

{$I procs6.pas} — { INTENSITY and UTILITY PROCEDURES } 


{PA} 


EEREKER AA ERE AS EE Se Fee eee eS KK KKK KKK KKK KKK KKEKAKKAKEKEKEKE 


* CALLED FROM: MAIN 

* CALLS TO: PUSH, POP, CALCULATE INTENSITY 

* DESC : Calculates the light intensity at a given intersection point 
* INPUT : A complete ray data record, the pointer to the light array 
. and the pointer to the object array 

* OUTPUT : The intensity at a given intersection point. If the input 
_ was a view ray then this intensity will be the intensity 
displayed. If the input ray is a reflected ray then this 
intensity is assigned to the I s field in the source ray. 

If the input ray 1s a refracted ray then this intensity is 
assigned to the | t field in the source ray. 


AMM K MM MMMM KK KK KK MK KOK KK KOR K Ke KOK oe oe KOK OR Ke OK KK KOK ROKK KKK KOK KKK KKK mk kh KKK KK 


* 
* 
* 


«x 


f 


procedure 12 (1 d__ : real; 
1 I tr,1 Iotepy lt tbh ©: real; 
1_[ sr, 1clisgy sb reas, 
i number of light sources : integer; 
i ambient r,i ambient g,1 ambient b_ :: real; 
i K ar,i K.ag,i K ab : real; 
1K g,i K Sg,1 K sb ; real; 
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i K tr,i K tg,i K tb : real; 
i K dr,i K dg,i K db : Teal: 
i ray type : raytype; 

1 ray top :ray_ ptr; 

1 light top: light ptr; 

i_obj ptr : object pir: 

1 0d) phong exp integer: 

1 ray vector x, 

1 ray vector y, 

1 ray vector z : real; 

1 intersection x, 

1 intersection y. 

1 intersection z : real; 

i surface normal x, 

i “surface “normal - ys 

1 “surface “normal 2 : real: 


var jo ray generation number é 
var o intensity red, 
oO intensity green, 
oO intensity _ blue : : reat); 


var 


integer; 


{ SET UP TEMPORARY AREA TO HOLD RAYS POPPED FROM STACK } 


templ ray type : raytype; 
templ ray origin x, 
templ ray origin y, 
templ ray origin z: 
templ ray vector x, 
templ ray vector y, 
templ ray vector z: real; 
templ ray stype =: raytype; 
templ intersection flag : boolean; 
templ _obj_ idx, 

templ subobj_idx, 

templ cpart idx, 


real; 


templ polygon idx : Integer; 
templ intersection x, 

temp] intersection y, 

templ | "intersection azn Teal, 


templ | ds: real; 

templ | J _tr, temp! I tg, templ I tb: 
templ Ae _sr, templ dT _sg, templ | ‘ “sb: 
temp2 ray type -: raytype: 
temp2 ray origin x, 
temp2 ray origin y, 
temp2 ray origin z : real; 
temp2 ray vector x, 
temp2 ray vector y, 
temp2 ray vector z : real; 

temp2 ray stype  :raytype; 
temp2 intersection flag : boolean; 


real; 
real: 
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temp2 obj idx, 

temp2 subobj idx, 

temp2 cpart idx, 

temp2 polygon idx : Integer; 

temp2 intersection x, 

temp2 intersection y, 

temp2 intersection z_ : real; 

temp2 | ds: real; 

temp2 _ J _tr, temp2 I tg, temp2 I tb: real; 
temp2_ alee _sr, temp2 | al _sg, temp2 | el “sb: real; 


begin 


{ CALCULATE INTENSITY OF THE RED COMPONENT OF LIGHT } 
calculate intensity( red, 
1d. 
i ambient r, 
1 (te, 
1 esr, 
1 _Tfay vector oe _Tay vector yl _Tay vector Z, 
1 _number _of light _sources, 
i K ar, 
i Kase 
1 -K tr, 
1 _K _dr, 
1 _obj _phong exp, 
1 intersection x, 
i intersection y, 
1 intersection 2, 
i surface _normal JX, 
1 "surface “normal - ny 
i surface “normal 4, 
i obj ptr, 
1 light top, 
o intensity red ); 


e 


{ CALCULATE INTENSITY OF THE GREEN COMPONENT OF LIGHT } 
calculate intensity( green, 
Leas 
1 ambient g, 
i I tg, 
1 Daege 
| ray vector X.1 ray vector y,1 ray vector 2, 
i number of light sources, - - ’ 
1 _K _ag, 
i K sg, 
i K tg, 
i K dg, 
1 obj phong exp, 
1 intersection x, 
1 intersection y, 
1 intersection 2, 
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1 surface normal x, 
i surface normal y, 
i surface norm al z, 
1 obj ptr, 

1 light top, 

o intensity green ); 


{ CALCULATE INTENSITY OF THE BLUE COMPONENT OF LIGHT } 
calculate intensity( blue, 
id, 
i ambient _b, 
oat. 
i I sb, 
1 ray vector x, 1 ray vector y, 1 ray vector 2, 
1 number of light sources, 
i K ab, 
1 K sb, 
i K tb, 
i K db. 
i obj phong exp, 
1 intersection x, 
1 intersection y, 
1 intersection 2, 
1 surface normal x, 
i surface normal y, 
i surface normal z, 
i obj ptr, 
i light top, 
o intensity blue ); 


if (i ray type = view) then begin 


THE VIEW RAY IS ALWAYS THE LAST RAY ON THE STACK WHEN IT IS POPPED THE 
INTENSITY DETERMINED FOR IT IS PASSED BACK INTO MAIN FOR DISPLAY. 


} 


{ nothing } 
end 
else begin 


if (i ray type = reflected) then begin 


sever | VistisMnm REFLECTED RAY THEN JUST ONE RAY NEEDS TO*BE POPPED TO 
GAIN ACCESS TO THE SOURCE RAY. 
j 

pop (templ ray type, 
templ ray origin x, templ ray origin y, templ ray origin z, 
templ ray vector x, templ ray vector y, templ ray vector z, 
templ ray stype, 
templ1 intersection flag, 
temp1 obj idx, 
templ subob) idx, 
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templ cpart idx, 

templ polygon idx, 

templ intersection x, temp] intersection y, templ intersection z, 
templ d, 7 
tempi I tr, templ_I tg, temp1 I tb, 

templ I sr, temp1 I sg, temp! I sb, 

ray top iG 


{ SET INTENSITY IN SOURCE RAY } 
templ I sr:=o intensity red; 
templ I sg :=o intensity green; 
temp1 I sb:=o intensity blue; 


{ RESTORE STACK } 

push( templ ray type, 
templ ray origin x, templ ray origin y, temp] ray origin z, 
templ ray vector x, templ ray vector y, templ ray vector Z, 
templ ray stype, 7 
temp! intersection flag, 
templ obj idx, 
temp1 subobj idx, 
templ cpart idx, 
templ polygon idx, 
templ intersection x, templ intersection y, templ intersection z, 
templ d, - 
templ I tr, templ I tg, temp1 I tb, 
temp1 I sr, temp1 I sg, temp1 I sb, 
ray top ); 


end 
else begin 


SINCE THIS IS THE REFRACTED RAY TWO RAYS MUST BE POPPED TO GAIN ACCESS 
TO THE SOURCE RAY. 
} 
pop (templ ray type, 

templ ray origin x, templ ray origin y, templ ray origin z, 
templ ray vector x, templ ray vector y, templ ray vector z, 
templ ray stype, 7 
templ intersection flag, 
templ obj idx, 
templ subobj idx. 
templ cpart idx, 
tempi! polygon idx, 
templ intersection x, tempi intersection y, tempi intersection 2, 
templ d, : 
templ I tr, templ I tg, temp1 I tb, 
templ I sr, templ I sg, temp! I sb, 
ray top IF 


pop (temp2 ray type, 
temp2 ray origin x, temp2 ray origin y, temp2 ray origin z, 


110 


temp2 ray vector x, temp2 ray vector y, temp2 ray vector 2, 
temp2 ray stype, 

temp2 intersection flag, 

temp2 obj idx, 

temp2 subob) idx, 

temp2 cpart idx. 

temp2 polygon idx, 

temp2 intersection xX, temp2 intersection y, temp2_ intersection 2, 
temp2 d, 

temp2 I tr, temp2 I tg, temp2 I tb, 

temp2 I sr, temp2 I sg, temp2 I sb, 

ray top}; 


{ SET INTENSITY IN THE SOURCE RAY. } 
temp2 I tr :=o intensity red; 
temp2 I tg:=o intensity green; 
temp2 I tb:=o intensity blue; 


{ RESTORE STACK. } 
push( temp2 ray type, 


temp2 ray origin x, temp2 ray origin y, temp2 ray origin 2, 
temp2 ray vector x, temp2 ray vector y, temp2 ray vector 2, 
temp2 ray stype, 

temp2 intersection flag, 

temp2 ob) idx, 

temp2_ subob) idx, 

temp2 cpart idx, 

temp2_ polygon idx, 

temp2 intersection x, temp2 intersection y, temp2 intersection 2, 
temp2 d, 

temp2 I tr, temp2 I tg, temp2_I tb, 

temp2 I sr, temp2 I sg, temp2 I sb, 

ray_top } 


push( templ ray type, 


end; 


templ ray origin x, templ ray origin y, templ ray origin z, 
templ ray vector x, templ ray vector y, templ ray vector 2, 
templ ray stype, 

templ intersection flag, 

templ obj idx, 

templ subob)j idx, 

templ cpart idx. 

temp] polygon ox, 

templ intersection <, tempi intersection y, temp! intersection 2z, 
templ d, 

templ I tr, templ I tg, templ I tb, 

templ I sr, templ I sg, templ I sb, 

ray top ); 


SINCE THE RAY SENT INTO THIS PROCEDURE IS NO LONGER NEEDED IT IS 
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DISCARDED, HENCE THE RAY COUNT NEEDS TO BE DECREMENTED. 
} 
io ray generation number := 10 ray generation number - 1; 
end; 
end; 


{ PA} 


KKK KKAKKHKKKKKKKKKKKHKKKKKAKAKEKKKEKE 16 KH KKKAKKKKHKKKKKKKKKKAKRKKKKKKKKKK KS 


* CALLED FROM : MAIN 

* CALLS TO : PUSH, POP, CALCULATE INTENSITY .. 
* DESC: Used to calculate intensity at node of stack storage is 

exceeded. 

* INPUT : Incomplete ray data record, this ray can not be continued 

. because there is no room on the stack for it. Also pointers 

to the light array and object array. 

* OUTPUT : Intensity, I t or I s, is set in source ray. 


2k RO RR i a aki ia ka ak ak ak aad akc ak ak ak: ae 


} 


* 


procedure I6 (i d : real; 
1 I tr, 1 ite ieee real; 
1 1 sr, i Iisewigeesb> 6: real: 
1 number of light sources : integer; 
i ambient r,i ambient g,1 ambient b : real; 
i K ar,i K ag,i K ab : real; 
me st; | Kasei sb : real; 
mK tr, 1 Kitg ko tb : real; 
i K dr, i Kodgiagk db sesreal: 
1 ray type : raytype; 
1 ray top : ray ptr; 
1 light top  : light ptr; 
1 obj ptr : object ptr; 
1 obj phong exp: integer; 
1 ray vector x, 
1 ray vector y, 
1 ray vector z : real; 
1 intersection x, 
1 intersection y, 
i intersection z_ : real; 
1 surface normal x, 
1 surface no rma! LY, 
i surface normal z : real: 
var o intensity red, 
o intensity green, 
o intensity blue : real ); 


var 

{ STORAGE FOR POPPED RAYS } 
templ ray type : raytype; 
templ ray origin x, 
templ ray origin y, 
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templ ray origin z: real; 
templ ray vector x, 

templ ray vector y, 

templ ray vector z: real; 

templ ray stype =: raytype; 

templ intersection flag : boolean: 

templ obj idx, 

templ subobj idx, 

templ cpart idx, 

templ polygon idx : integer; 

temp] intersection x, 

templ intersection y, 

templ intersection z : real; 

templ d_ : real; 

templ J _tr, templ I tg, templ I tb: real; 
temp] | Ts _sr, templ | “J _sg, templ | il sb: real; 


temp2 ray type : raytype; 

temp2 ray origin x, 

temp2 ray origin y, 

temp2 ray origin z: real; 
temp2 ray vector x, 

temp2 ray vector y, 

temp2 ray vector 2z: real; 
temp2 ray stype  : raytype; 

temp2 intersection flag : boolean; 

temp2_ ob) idx, 

temp2 _subobj - idx, 

temp2 cpart idx, 

temp2 polygon idx : integer; 

temp2 intersection x, 

temp2 intersection y, 

temp2 intersection 2 : real; 

temp2d» : real; 

temp2_ I _tr, temp2 I tg, temp2 I tb: real; 
temp2 | lw: _sr, temp2__ a _sg, temp2 | a “sb : real; 


begin 
{*** COULD INSERT TREE EXTENSION PROCEDURE ***} 


{ CALCULATE RED COMPONENT OF LIGHT } 
calculate intensity( red. 
id, 
1 ambient r, 
1 1 te 
1 hast; 
1 ray vector x,1 ray vector y,1i ray vector z, 
i number of light sources, 7 - 
1 Keir, 
i K sr, 
i K tr, 
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1 K dr, 

i obj phong exp, 

1 intersection x, 

1 intersection y, 

1 intersection 2, 

1 surface normal x, 
1 surface normal y, 
i surface normal z, 
1 obj ptr, 

i light top, 

o intensity red ); 


{ CALCULATE GREEN COMPONENT OF LIGHT } 
calculate intensity( green, 
i d, 
1 ambient g, 
i I tg, 
i I sg, 
1 ray vector x,1 ray vector y,1 ray vector z. 
1 “number _of light sources, 
i K ag, 
i K sg, 
1 K tg, 
1_K dg, 
1 obj phong exp, 
1 intersection x, 
1 intersection y, 
1 intersection z, 
1 _surface Pea ox 
1 _surface “normal | BY 
1 "surface “normal nz; 
i obj ptr, 
1 light top, 
o intensity green); 


{ CALCULATE BLUE COMPONENT OF LIGHT } 
calculate intensity( blue, 
i d, 
1 ambient b, 
i I tb, 
1 I sb, 
1 _Tay vector _X, 1 _Tay vector y, 1 Tay vector 2, 
1 number of light sources, 
i K ab. 
1 K_ sb, 
i K tb, 
i K db, 
1 obj phong exp, 
i intersection x, 
1 intersection y, 
1 intersection z, 
1 surface normal x, 
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\_ surface normal y, 
1 surface normal 2, 
1 obj ptr, 

1 light top, 

o intensity blue); 


if (i ray type = reflected) then begin 


IF INPUT RAY IS A REFLECTED RAY THEN ONE RAY MUST BE POPPED TO GAIN 
ACCESS TO IT’S SOURCE RAY. 
} 

pop (temp1 ray type, 
tempi ray origin x, templ ray origin y, templ ray origin z, 
templ ray vector x, templ ray vector y, templ ray vector 2, 
templ ray stype, 
templ intersection flag, 
temp1 obj _idx, 
templ subob) idx. 
templ cpart idx, 
templ polygon idx, 
templ intersection x, templ intersection y, templ intersection 2, 
templ d. 
templ I tr, tempi I tg, tempi I tb, 
templ I sr, tempi I sg, temp1 I sb, 
1 ray_top 1 


{ SET I s INSOURCE RAY } 
templ I sr :=o intensity red; 
templ I sg :=o intensity green; 
templ I sb :=o intensity blue; 


{ RETORE STACK. } 

push( temp1 ray type, 
templ ray origin x, templ ray origin y, templ ray origin z, 
templ ray vector x, templ ray vector y, templ ray vector 2, 
templ ray _ stype, 
templ intersection flag, 
templ obj idx, 
templ subob) idx, 
templ cpart idx, 
templ polygon idx. 
temp] “intersection _X, tempi intersection y, templ intersection 2, 
tempi d. 
templ I| tr, templ I tg, tempi I tb, 
templ I sr, templ I sg, temp1 I sb, 
ray top ); 


end 


else begin 


{ 
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[F [INPUT RAY IS A REFRACTED RAY THEN TWO RAYS MUST BE POPPED FROM THE 
STACK TO GAIN ACCESS TO IT’S SOURCE RAY. 


} 

pop (templ ray type, 
templ ray origin x, templ ray origin y, templ ray origin z, 
templ ray vector x, templ ray vector y, templ ray vector 2, 
temol ray stype, 
templ intersection flag, 
templ obj idx, 
templ subobj idx, 
templ cpart idx, 
templ polygon idx, 
templ intersection x, temp1 intersection y, temp1 intersection z, 
templ d, 
templ I tr, templ I tg, temp1 I tb, 
templ I sr, temp1 I sg, templ I sb, 
ray top _—); 


pop {temp2 ray type, 
temp2 ray origin x, temp2 ray origin y, temp2 ray origin z, 
temp2 ray vector x, temp2 ray vector y, temp2 ray vector 2, 
temp2 ray stype, 
temp2 intersection flag, 
temp2 obj idx, 
temp2 subobj idx, 
temp2 cpart idx, 
temp2 polygon idx, 
temp2 intersection x, temp2 intersection _y, temp2_ intersection 2, 
temp2 d, 
temp2 I tr, temp2 I tg, temp2 I tb, 
temp2 I sr, temp2 I sg, temp2 I sb, 
ray top): 


{ SET THE I t FIELD IN THE SOURCE RAY } 
temp2 I tr:=o intensity red; 
temp2 I tg :=o intensity green; 
temp2 I tb :=o intensity blue; 


{ RESTORE THE STACK } 

push( temp2 ray type, 
temp2 ray origin x, temp2 ray origin y, temp2 ray origin 2, 
temp2 ray vector x, temp2 ray vector y, temp2 ray vector 2, 
remp2 ray stype, 
temp2 intersection ilag, 
temp2 obj idx, 
temp2 subobj idx, 
temp2 cpart idx, 
temp2 polygon idx, 
temp2 intersection x, temp2 intersection y, temp2 intersection z, 
temp2 d, 
temp2 I tr, temp2 I tg, temp2 I tb, 
temp2 I sr, temp2 I sg, temp2 I sb, 
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ray top ; 


push( temp1 ray type, 
templ ray origin x, templ ray origin y, templ ray origin z, 
templ ray origin x, templ ray origin y, templ ray origin z, 
templ ray stype, 
temp! intersection fag, 
templ obj idx, 
templ_ subobj idx, 
templ_ cpart idx, 
temp1 polygon idx, 
tempi intersection x, templ intersection y, temp! intersection 2, 
templ d, 
templ I tr, templ I tg, temp1 I tb, 
templ I sr, templ I sg, temp1 I sb, 
ray top i 

end; 


end: 


{PA} 


{ EMER EERE EKKEKEEREEE SEEK EKREK SEEKER EKER KEE KEK RE EK EEK EK EEE KEK KKK KE 

REERERKREKEKTKRKKEKSEEKKERKEKEKEKKEEKRKEEEKSE MAIN KEKKEKEKKKRKKEKEKEKEKEKKREKKEE KEE KK KK 
{ KEKKEKEKKEEKEEKK EKER KEK KKK KEK KKK EK KEKE KKK KKK EEK KK KEK KK KEKE KKK EK KEKE KES \ 
begin 

assign (sysin, picture5.pas’); 

reset (sysin); 


{ OUTPUT FILE } 
assign (outfile, pic0.dta’); 
rewrite (outfile); 


{ KKK KEKE KEKE KEK KEKE KEKE KEKE EEE EEE EEK KEKE EEE KEKE KEKE KE KEKE KEKE KEE KE KES } 


{ SET UP COUNTERS } 
obj) cntr) :=1; 
subobj cntr := 1; 
hight cntr := 1; 


cpart cntr := 1; 
poly cntr  -:=#b; 
verlice  Cnlr = — 15 


{ CREATE POINTERS TO RECORDS } 
new(obj curr); 
new(subob) curr); 
new(cpart_ curr); 
new(light current); 
new(poly curr); 
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{ *** READ IN DATA FILE *** } 


{ READ IN PICTURE RECORD } 
readin (sysin, picture.view position x); 
readin (sysin, picture.view position y); 
readIn (sysin, picture.view position 2); 
readIn (sysin, picture.background coior +}; 
readin (sysin, picture.background | “color _g); 
readIn (sysin, picture.background | “color ab), 
readIn (sysin, picture. screen max x); 
readIn (sysin, picture.screen_max_y); 
readln (sysin, picture.ambient_r); 
read]n (sysin, picture.ambient_g); 
readIn (sysin, picture.ambient_b); 
readIn (sysin, picture.no_ zero); 
read]n (sysin, picture.global refraction index); 
readln (sysin, picture.num _lights); 
picture.lights := light current; 


while picture.num_ lights > 0 do begin 


{ READ IN LIGHT DATA } 

readin (sysin,light current “[light _cntr].1 _); 
readIn (sysin,light current “(light cntrj.I_g); 
readIn (sysin,light current “light saith mb); 
readIn (sysin,light current “[light cntr].light x); 
readIn (sysin,light current “[light cntr|.light_y); 
readIn (sysin,light current “[light cntr].light_z); 
readln (sysin, light current” [light _cntr].dimension1); 
readIn (sysin,light_ current “[light cntr].dimension2); 
light _cntr := hight cntr + 1; 
picture.num lights := picture.num lights - 1; 

end; 

read]n (sysin, picture.num_objs); 

object loop cnt := picture.num_objs; 

picture.objects := obj curr; 


while object loop cnt > Odo begin 

{ READ IN OBJECT DATA } 
readIn (sysin, obj curr “{obj_cntr].opcode); 
readin (sysin, obj curr“!obj cntr|.obj bsphere radius); 
readin (sysin, ob) eee _cntr}.obj _bsphere | et 
readin (sysin, obj curr ~/ob) cntrj.obj) bsphere _y); 
readin (sysin, obj curr “|obj cntrj.ob) _bsphere _2); 
readin (sysin, obj curr“[obj cntr].num_ sub objects); 
subobj loop cnt := obj curr “[obj _entr]. num sub objects; 
picture. objects {obj _entr]. sub objects := subobj curr; 


while subobj loop cnt > 0do begin 
{ READ IN SUBOBJECT DATA } 
readIn (sysin, subobj curr “|subobj cntr].subobj type); 
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readin (sysin, subob) curr *'subod}j _cntr|.sub bsphere radius); 

readIn (sysin, subobj curr” subobj cntr|.sub bsphere x); 

readIn (sysin, subobj curr ~'subobj _entr). sub bsphere y); 

readIn (sysin, subobd) curr jsubobj cntrj.sub_bsphere z); 

readIn (sysin, subobj curr”/subobj cntr].num_common_parts); 

cpart loop cnt := su bob} Scunr *{subob} _entr]. num common parts; 
dicture.objects ~ ‘obj cntr.. subd objects” isubobj centri. .common parts := 
cpart curr; 


while cpart loop cnt > Odo begin 
{ READ IN COMMON PART DATA } 
readIn(sysin, cpart curr” cpart cntr!.K ar); 
readln(sysin, cpart curr ~|cpart _entr]. K_ag); 
readIn(sysin, cpart curr“!cpart cntr].K ab); 
readIn(sysin, cpart curr “lepart _entr].K dr); 
readIn(sysin, cpart curr“{cpart cntr].K dg); 
readIn(sysin, cpart_curr “[cpart _entr|.K_ db); 
readIn(sysin, cpart curr ~ cpart _entr]. K sr); 
readInjsysin, cpart curr” cpart cntr}.K sg); 
readIn(sysin, cpart curr ,cpart cntr|.K sb); 
readIn(sysin, cpart_curr”|cpart cntr].K tr); 
readIn(sysin, cpart curr”|cpart_cntr|.K_ tg); 
readIn(sysin, cpart curr” cpart__ entr!.K tb); 
readIn(sysin, cpart curr” cpart cntrj.obj refraction index); 
readIn(sysin, cpart curr” cpart_cntr].obj _phong _exp); 


{ CHECK TO SEE IF SUBOBJECT IS A SPHERE OR A POLYGONAL OBJECT } 
if (subobj_curr*{subobj cntr].subobj type = 1) then begin 


readIn(sysin, cpart curr “[cpart cntr].num_polygons); 
poly _loop | cnt := cpart curr “lcpart _entr]. num _polygons; 
picture.objects ~ [obj _entr]. sub objects *[subobj _entr}. 
common parts *|cpart cntr].polygons := poly curr; 


while poly loop cnt > Odo begin 
{ READ IN POLYGON DATA } 
readIn(sysin, poly curr*[poly cntr].num_vertices); 
vertice loop cnt := poly curr*|poly cntr].num_vertices; 


while vertice loop cnt > 0 do begin 
{ READ IN VERTICE DATA } 

readIn(sysin. poly curr“!poly cntr].vertice x/vertice cntr!); 

readInisysin, poiy curr “!poiy cntrj.vertice y|vertice “entr}); 

readln(sysin. poly curr” poly _cntr}.vertice _aivertice “cntr)): 
vertice cntr := vertice cntr + 1; . 
vertice loop cnt := vertice loop cnt - 1; 

end; 

vertice cntr := 1; 

readIn(sysin, poly curr *|poly cntr|.surface normal x); 

read]n(sysin, poly curr” [poly _cntr]. surface normal _y); 

readln(sysin, poly curr*{poly cntr|.surface normal z); 

poly cntr := poly cntr + 1; 
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poly loop cnt := poly loop cnt - 1 


end; 
poly cntr := 1; 


new(poly curr); 


end; { [F } 
cpart cntr :=cpart cntr + 1; 
cpart loop cnt := cpart loop cnt - 1 


end; 

cpart cntr:= 1; 

new(cpart curr); 

subobj cntr := subob) cntr + 1; 

subobj loop cnt := subobj loop cnt - 1 


end; 
new(subobj curr); 
subobj cntr := 1; 


obj cntr := obj cntr + 1; 
object loop cnt := object loop cnt - 1 


end; 


RRR RK KKK KEKE REE KR KEK KEK KK KR KEK KKK KKK KKK KKK KKK KKK KKK EK KEK KEK KEK KKK KKK KEK EK \ 


{ SET RAY STACK POINTER } 


ray top) := nil; 
pixel z := initial pixel z; 


{ RASTER SCAN LOOP } 


for pixel _y := 1 to picture.screen_max_y do begin 
for pixel _x := 1 to picture.screen max _ x do begin 


{ DETERMINE VIEW RAY DIRECTION } 
x := pixel x - picture.view position x; 
y := pixel y - picture.view position _y; 
z:= pixel z- picture.view position 2; 
dist := sart{sar(x) + sar{y) + sqr(z)); 
{ CONVERT I[T TO A UNIT VECTOR } 
unIcx = x fedise 
unity := y / dist; 
unitz := z / dist; 


{ INITIALIZE VIEW RAY } 
initial ray type := view; 
initial ray origin x := pixel x; 
initial ray origin y := pixel y; 
initial ray origin z:= pixel 2; 
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initial ray vector x := unitx; 
initial ray vector y := unity; 
initial ray vector z := unitz; 
initial ray stype := none; 
initial intersection ine := false; 
initial “obj | idx s= 0; 
initial subobj idx at, 
initial _cpart idx a0). 


initial polygon idx s= 0; 


initial intersection x =a.) 
initial | _intersection y = := 0.0; 
initial » _intersection 2 — 00: 
initial d  := 0.0; 
initial I tr —— 0.0; 
initial I tg := 0.0; 
initial I tb := 0.0; 
initial I sr := 0.0; 
initial I sg := 0.0; 


inivaas 1 sb := 0.0: 
ray generation number := 0; 


push( initial ray type, 


initial ray origin x,initial ray origin y,initial ray origin 2, 
initial ray vector x,initial ray vector y,initial ray vector z, 


initial ray stype, 
initial intersection flag, 
initial obj idx, 

initial subob)j idx, 
initial cpart idx, 
initial polygon idx, 
initial intersection x, 
initial | intersection _y, 
initial | _intersection _2, 
initial Pal 


initial I _tr, initial I tg, initial I tb, 


initial LT: _sr, initial u _ sg, initial al “sb, 
ray top 1 


{ BEGIN RAY TRACING LOOP } 


repeat 


pop (current ray type. 


current ray origin x,current ray origin y,current ray origin 2, 
current ray vector x,current ray vector y,current ray vector 2, 


current ray stype, 
current intersection flag, 
current obj idx, 

current subobj idx, 
current cpart idx, 
current polygon idx, 
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current intersection x, 

current intersection y, 

current intersection 2, 

current d, 

current I tr, current I tg, current I tb, 
current I sr, current I sg, current I sb, 


ray top A 
if (current intersection flag) then begin 


IF THIS FLAG IS SET THEN THIS RAY HAS ALREADY BEEN THROUGH THE RAY 
TRACING PROCESS AND HAS HIT AN OBJECT. 
} 
cpart path := picture. 
objects “[current obj idx]. 
sub objects *[current subobj_ idx]. 
common Parts; 


subob) path := picture. 
objects *|current obj idx]. 
sub_ objects; 


if subob) path “[eurrent subobj idx|.subobj type = 0 then begin 


IF THE SUBOBJECT TYPE IS A SPHERE THEN THE SURFACE NORMAL AT THE POINT 
OF INTERSECTION MUST BE CALCULATED SINCE IT CAN NOT BE STORED. THE SURFACE 
NORMAL IS DETERMINED FOR THE PLANE TANGENT TO THE SPHERE AT THE 
INTERSECTION POINT. 
} 
surface normal x := current intersection x - 
subob} _path * [current | subobj _idx)]. 
sub bsphere x; 


surface _normal _y := current intersection y= 
subobj path * {current _subobj _idx]. 
sub bsphere y; 


surface normal z := current intersection z - 
subobj path “|current subobj idx]. 
sub bsphere z; 


dist := ( sqrt({ sqr(surface normal x) + 
sqr(surface normal y) + 
sqr(surtace normal 2)})); 


{ THIS RAY IS THEN CONVERTED INTO A UNIT VECTOR } 
surface normal x := surface normal x / dist; 
surface normal y := surface normal — _y / dist; 
surface “normal _2 := surface “normal _2/ dist; 
end 
else begin 
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IF THE SUBOBJECT IS A POLYGONAL OBJECT THEN THE SURFACE NORMALS FOR EACH 
OF THE POLYGONS OF WHICH IT IS COMPOSED IS RETRIEVED FROM IT’S RECORD. 
} 
surface normal x := cpart path “jcurrent cpart idx). 
polygons “{current polygon idx). 
surface normal x: 


surface normal y := cpart path “{current cpart idx}. 
polygons “(current polygon idx). 
surface normal y; 


surface normal z:=cpart path “[current cpart idx]. 
polygons ~ [current polygon idx]. 
surface normal 2; 
end; 


PROCEDURE FOR DETERMINING THE INTENSITY OF LIGHT AT EACH INTERSECTION 
POINT. 
} 
I2 (current d, 
current wil _tr, current I tg, current I tb, 
current iz _sr, current wl _sg, current za tb, 
picture.num lights, 
picture.ambient r, picture.ambient g, picture.ambient b, 
cpart path “[current | cpart idx].K ar, 
cpart path “(current cpart idx]. K _ag, 
cpart path “(current cpart_idx!.K ‘ab, 
cpart path “(current cpart idx|.K _dr, 
cpart path” [current _cpart _idx|. K _dg, 
cpart path “|current cpart idx].K _db, 
cpart path ~ [current cpart idx]. K sr, 
cpart path “|current cpart idx].K sg, 
cpart path “(current cpart idx].K sb, 
cpart path “|current cpart idx]. K _tr, 
cpart path “(current cpart idx]. K _tg, 
cpart path “(current cpart_idx].K tb, 
current ray type, 
ray top, 
picture.lights, 
picture.objects, 
cpart path ‘current cpart idx|.obj phong exp, 
current ray vector x. J 7 7 
current ray _vector wy, 
current ray vector 2, 
current intersection x, 
current intersection y, 
current _Intersection 4, 
surface _normal _X, surface normal y, surface normal z, 
eect len seumber, 
intensity red, 
intensity green, 
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intensitv blue}; 


end 
else begin 


{ CHECK FOR POSSIBLE INTERSECTIONS OF CURRENT RAY WITH OBJECTS IN SCENE } 
check for intersection(current ray vector x, 
- current ray vector y, ; 
current ray vector 2Z, 
current ray origin x, 
current ray origin y, : 
current ray origin 2, 
picture.objects, 
current intersection x, 
current intersection y, 
current intersection 2, 
current obj idx, 
current subobj idx, 
current cpart idx, 
current polygon idx, 
current intersection flag); 


{ SET UP PATHNAMES TO USE AS SHORTHAND } 
subobj path := picture. 
objects *{current obj idx]. 
sub objects; 


cpart path := picture. 
objects *{current obj idx]. 
sub objects {current subobj idx]. 
common parts; 


if subobj path “[current subobj idx].subobj type = 0 then begin 
surface normal x := current intersection x - 
subobj path “|current subobj idx]. 
sub bsphere x; 
surface normal y := current intersection y - 
subobj path “|current subobj_ idx]. 
sub bsphere y; 
surface normal z := current intersection z - 
subobj path “|current subobj idx]. 
sub bsphere 2: 
dist := ( sart( sqr(surface normal x) — 
sqr(surtace normal y) + 
sqr(surface normal z))); 


surface normal x := surface normal x / dist; 
surface normal y := surface normal y / dist; 
surface normal z := surface normal z / dist; 
end 
else begin 


124 


surface normal x :=cpart path “(current cpart idx]. 
polygons “{current polygon idx]. 
surface normal x; 


surface normal y :=cpart path “{current cpart idx]. 
polygons “[current polygon idxl. 
surface normal y; 


surface normal z :=cpart_path“|current cpart idx]. 
polygons “(current polygon idx]. 
surface normal 2; 
end; 


if (current intersection flag) then begin 


{ IF THERE HAS BEEN AN INTERSECTION THEN CONTINUE TRACING THE RAY } 


if (stack exceeded( ray generation number, 
maximum size of _stack)) then begin 


{ IF THE STACK IS ALREADY FULL THEN CALCULATE INTENSITY AT LAST NODE } 


I6 (current d, 
current J _tr, current I tg, current I tb, 
current cis _3r, current zi _sg, current a pcb 
picture. num lights, 
picture.ambient r, 
picture.ambient g, 
picture.ambient_b, 
cpart path “(current cpart idx].K ar, 
cpart path “|current cpart idx].K ag, 
cpart path “|current cpart idx|.K ab, 
cpart path “|current cpart_idx|.K dr, 
cpart path “|current cpart_idx|.K _dg, 
cpart path “|current_cpart_idx|.K db, 
cpart path “|current cpart idx). K : _sr, 
cpart path” [current cpart idx| K sg, 
cpart path “|current cpart idx].K_ sb, 
cpart path “(current cpart idx|.K tr, 
cpart path “{current cpart idx].K tg, 
cpart path “|current cpart idx|.K tb, 
current ray type, 
ray top, 
picture. lights. 
picture.opjects. 
cpart path current cpart idxi.ob} phong exp. 
current ray vector x, - 
current ray vector y, 
current ray vector 2, 
current intersection x, 
current intersection y, 
current intersection 2, 
surface normal x, 
surface normal y, 
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surface normal z, 
intensity red, 
intensity green, 
intensity blue) 
end 
else begin 


IF THERE WAS AN INTERSECTION AND THE STACK WAS NOT FULL THEN CALCULATE 
THE DISTANCE BETWEEN THE RAY’S ORIGIN AND POINT OF INTERSECTION AND PLACE 
THE RAY BACK ON THE STACK. 


} 


current d = (sqrt(sqr(current intersection x - 
current ray origin x)) + 
(sqr(current intersection x - 
current ray origin x)) + 
(sqr(current intersection x - 
current ray origin x))); 


push( current ray type, 
current ray origin x, 
current ray origin y, 
current ray origin 2, 
current ray vector x, 
current ray vector y, 
current ray vector 2, 
current ray stype, 
current intersection flag, 
current obj idx, 
current subob) idx, 
current cpart idx, 
current polygon idx, 
current intersection x, 
current intersection y, 
current intersection 2, 
current d, 
current I tr, 
current | tg, 
current I tb, 
current I sr, 
current I tg, 
current { tb, 
ray top ); 


{ DETERMINE IF A REFLECTED RAY WAS CREATED AND IF SO CALCULATE IT } 
calculate reflected ray (current ray vector x, 

current ray vector y, 

current ray vector Z, 

surface normal x, 

surface normal y, 

surface normal z, 

reflected ray x, 
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reflected ray y, 
reflected ray 2, 
reflected ray +); 


{ DETERMINE IF A REFRACTED RAY WAS CREATED AND IF SO CALCULATE IT } 


calculate retracted ray (current ray vector x, 
- current ray vector y, 

current ray — “vector _%, 
surface “normal ox, 
surface “normal - EY: 
surface “normal z 
cpart path” reuntertt _cpart idx]. 
obj refraction index, 
picture.global | refraction _index, 
refracted ray_x, 
refracted ray y, 
refracted ray 2, 
refracted ray +); 


ray generation number := ray generation number + 1; 


if reflected ray then begin 


IF A REFLECTED RAY WAS CREATED THEN INITIALIZE IT AND PUSH IT ON THE 


STACK. 
} 


source ray type := current ray type; 


dist := ( sqrt(sqr(reflected ray x) + 
sqr(reflected ray y) + 
sqr(reflected ray _2))); 


{ CONVERT REFLECTED RAY TO A UNIT VECTOR } 
unitx := reflected ray x / dist; 
unity := reflected ray y / dist; 
unitz := reflected ray z / dist; 


initial ray type := reflected; 
initial ray origin x := current intersection x; 
imitial ray origin y := current intersection y; 
initial ray origin 7 := current intersection 2: 
initial ray vector X< := uMmitx; 
initial ray vector y := unity; 
initial ray vector z := unitz; 


initial ray stype := source ray type; 
initial intersection flag := false; 
initial obj idx a) 

initial _subobj | idx := 0; 

initial cpart | idx := 0; 

initial polygon idx = 0; 
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initial intersection x = := 0.0; 
initial intersection y := 0.0; 
initial intersection 2 := 0.0; 
initial d a 
Initial ltr ee — Oee 


initial I te  := 0.0: 
initial ( tb = 0.0; 
initial [ sr = 0.0; 
initial [ sg  := 0.0; 
initial I sb = := 0.0; 


push( initial ray type, 
initial ray origin x, 
initial ray origin y, 
initial ray origin 2, 
initial ray vector x, 
initial ray vector _y, 
initial ray vector 2z, 
initial ray stype, 
initial intersection flag, 
initial ob) idx, 
initial subobj idx, 
initial cpart idx, 
initial _polygon_idx, 
initial intersection x, 
initial intersection y, 
initial intersection 2, 
initial d, 
initial I tr, initial I tg, initial I tb, 
initial I sr, initial I sg, initial I sb, 
ray top) 


end; 


if refracted ray then begin 


IF A REFRACTED RAY WAS CREATED THEN INITIALIZE IT AND PUSH IT ON THE 
STACK. 


source ray type := current ray type; 


dist := ( sqrt(sqr(refracted ray x} + 
sqr(refracted ray y) + 
sqr(refracted ray z))); 


{ CONVERT IT TO A UNIT VECTOR } 
unitx := refracted ray x / dist; 
unity := refracted ray _y / dist; 
unitz := refracted ray z / dist; 
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initial ray type 
initial ray origin x: 
initial ray origin y : 
initial ray origin z: 


:= refracted; 


current intersection x; 
current intersection y; 


current intersection 2; 


initial ray vector x := unitx; 
initial ray vector y := unity; 
imitial ray vector z := unitz; 
initial ray stype := source ray type; 
initial intersection flag := false; 
initial obj idx 2= 70>: 
initial subob) idx so) 
initial cpart idx 10): 
initial polygon idx =10; 
initial intersection x ai: 
initial intersection y  := 0.0; 
initia] intersection z =—oOs0. 


initial d = 0.0; 
initial I tr := 0.0; 
initial I tg wrt). 
initial I tb == 30; 
initial I sr = i070: 
initial I sg = ).0: 
initial I sb = =:= 0.0; 


push( initial ray type, 
initial ray origin x, 
initial ray origin y, 
initial ray origin z, 
initial ray vector x, 
initial ray vector y, 
initial ray vector z, 
initial ray stype, 
initial intersection flag, 
initial _obj idx, 
initial subobj idx, 
initial cpart idx, 
initial _polygon_ idx, 
initial intersection x, 
initial intersection _y, 
initial intersection z, 
initial d, 
initial ci _tr, initial I tg, initial I tb, 
initial LT: _sP. initial I _38, initial a eb. 
ray top 


}e 
e 
} 4 


end; 


end; 
end 
else begin 
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if (current ray type = view) then begin 


IF THERE WAS NO INTERSECTION AND THE CURRENT RAY IS THE VIEW RAY THEN SET 
THE OUTPUT INTENSITY TO THE BACKGROUND INTENSITY. 


intensity red := picture.background color r; 
intensity green := picture.background color g; 
intensity blue := picture.background color _b; 
end; 
end; 
end; 


until (stack empty(ray_top)); 
{ OUTPUT THE FINAL INTENSITY } 
display pixel( intensity red, 
intensity green, 
intensity blue, 


pixel x, 
pixel y, 
pixel z ); 
end 
end; 


close (sysin); 
close (output); 


end. { MAIN } 
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APPENDIX B — INPUT FILE 


100 x, view position { PICTURE } 
100 /y/ 

1000 /2/ 

0.0 /red/ background light 

0.0 /green/ 

1.0 /blue/ 

200 /x/ screen 

200 /y/ 

1.0 /red/ ambient intensity 

1.0 /green/ 

1.0 /blue/ 

1.0 /global refraction index/ 

1 num lights/ 

1.0 ,red/ intensity of light source 

10 /green/ 

1.0 /blue/ 

0.0 /x/ position of light source . 
20.0 ;y/ 

0.0 /2/ 


0.0 /dimension1/ 

0.0 /dimension2/ 

3  /num_ objects/ 

9999 /opcode/ { OBJECT 1} 
35.0 /radius of object’s bounding sphere/ 
110.0 /x/ center of bounding sphere 

0.0 /y/ 

-50.0 /z/ 

1  /number of subobjects/ 

1 /subobject type/ {SUBOBJECT 1} 
35.0 /radius of subobjects bounding sphere/ 
110.0 /x/ center of bounding sphere 

0.0 /y/ 

-50.0 /z/ 

1  /num-common-parts/ { COMMON PART 1 } 
0.8 /Ka-red/ ambient coefficient 

0.0 /Ka-green/ 

0.0 Ka-blue, 

0.3 Kd-red/ diffuse coefficient 

0.0 ,Kd-green;/ 

0.0 /Kd-blue/ 

0.8 /Ks-red/ specular coefficient 

0.8 /Ks-green/ 

0.8 /Ks-blue/ 

0.0 /Kt-red/ transmission coefficient 
0.0 /Kt-green/ 

0.0 /Kt-blue/ 
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0.0 /obj-refraction-index/ 
200 /obj-phong-specular-exponent / 
6 /num-polygons/ 

4 /num-vertices/ { POLYGON 1 } 
90.0 /polyl ptl/ 

20.0 

-30.0 

90.0 /poly! pt2/ 

-20.0 

-30.0 

130.0 /polyl pt3/ 

-20.0 

-30.0 

130.0 /polyl pt4/ 

20:0 

-30.0 

0.0 /poly1 surface normal/ 
0.0 

FO 

4  /num-vertices/ { POLYOGN 2 } 
90.0 /poly2 pt1/ 

-20.0 

-20.0 

90.0 /poly2 pt2/ 

-20.0 

-70.0 

130.0 /poly2 pt3/ 

-20.0 

-70.0 

130.0 /poly2 pt4/ 

-20.0 

-70.0 

0.0 /poly2 surface normal/ 
-1.0 

0.0 

4  /num-vertices/ {POLYGON 3} 
130.0 /poly3 pt1/ 

20.0 

-70.0 

130.0 /poly3 pt2/ 

-20.0 

-70.0 

90.0 Mpoly2 pts/ 

-20.0 

-70.0 

90.0 /poly3 pt4/ 

20.0 

-70.0 

0.0 /poly3 surface normal/ 
0.0 

-1.0 

4 /num-vertices/ { POLYGON 4 } 
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130.0 /poly4 pt1/ 

20.0 

-30.0 

130.0 /poly4 pt2/ 

20.0 

-70.0 

90.0 ,poly4 pt3/ 

20.0 

-70.0 

90.0 /poly4 pt4/ 

20.0 

-30.0 

0.0 /poly4 surface normal/ 
1.0 

0.0 

4 /num-vertices/ { POLYGON 5 } 
130.0 /poly5 pt1/ 

20.0 

-30.0 

130.0 /polyd pt2/ 

-20.0 

-30.0 

130.0 /polyd pt3/ 

-20.0 

-70.0 

130.0 /poly5 pt4/ 

20.0 

-70.0 

10 /poly5 surface normal/ 
0.0 

0.0 

4 /num-vertices/ { POLYGON 6 } 
90.0 /poly6 pt1/ 

20.0 

-70.0 

90.0 /poly6 pt2/ 

-20.0 

-70.0 

90.0 /poly6 pt3/ 

-20.0 

-30.0 

90.0 /poly6 pt4/ 

20.0 

-20.0 

-1.0 /poly6 surface normal/ 
0.0 

0.0 

9999 /opcode/  Obieer 2 >** 
175 /radius of the objects bounding sphere/ 
100.0 /x/ center point of the bounding sphere 
-100.0 /y/ 

-100.0 /z/ 
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1 /num-subobjects/ 

1 /subobject-type/ {SUBOBJECT 1 } 

175 _—/radius of subobjects bounding sphere/ 

100.0 /x/ center of bounding sphere 

-100.0 /y/ 

-100.0 /z/ 

1 /number common parts/ { COMMON PART 1 } 
0.0 /red/ Ka ambient coefficient 


0.7 /green/ 

0.0 /blue/ 

0.0 /red/ Kd diffuse coefficient 
0.7 /green/ 

0.0 /blue 

0.8 /red/ Ks specular coefficient 
0.8 /green/ 

0.8 /blue/ 

0.0 /red/ Kt transmission coefficient 
0.0 /green/ 

0.0 /blue/ 


0.0 /objects refraction index/ 

200 /Phong’s specular exponent/ 

1 /number of polygons/ 

4 /number of vertices/ { POLYGON 1 } 
0.0 /polyl ptl/ 

20.0 

-200.0 

0.0 /poly1 pt2/ 

0.0 

0.0 

200.0 /poly1 pt3/ 

0.0 

0.0 

200.0 /polyl pt4/ 

20.0 

-200.0 

0.0 /polyl surface normal/ 

0.99 

0.1 

9999 /opcode/ at OBIEC T 3 ** se 
40 /radius of objects bounding sphere/ 
140 /x/ center of bounding sphere 


30—Ss /y/ 
-150 /za; 
1 /number of subobjects/ 


0 /subobject type/ 

40 /radius of subobjects bounding sphere/ 

140 /x/ center of bounding sphere 

30 /y/ 

-150 /z/ 

1 /number of common parts/ { COMMON PART 1 } 
0.5 /red/ Ka ambient coefficient 

0.0 /green/ 


134 


0.5 
0.5 
0.0 
0.5 
0.8 
0.8 
0.8 
0.0 
0.0 
0.0 
0.0 
200 


/blue/ 

/red/ Kd_ diffuse coefficient 
/green/ 

/blue/ 

/red/ Ks specular coefficient 
/green/ 

blue/ 

/red/ Kt transmission coefficient 


/green/ 
/blue/ 


/refraction index for object / 
/Phong’s specular exponent / 
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